ESTABELECIMENTO E INTERPRETAÇÃO FÍSICA DA TRANS- 
FORMAÇÃO DE LORENTZ. 


Physical establishment and interpretation of the Lorentz trans- 
formation. 


M. ABREU FARO 


INTRODUÇÃO AOS CONJUNTOS VAGOS. 
Introduction to fuzzy sets. 
ANTÓNIO GOUVÊA PORTELA 
EDUCATION AND RESEARCH IN MECHANICS IN THE FE- 
DERAL KREPUBLIC OF GERMANY. 
Ensino e investigação em mecânica na República Federal Alemã. 
ERNST BECKER 
CONTINUOUS HYDROLOGICAL TIME-SERIES DISCRETI- 
ZATION. | 
Discretização de séries hidrológicas contínuas. 
LUÍS VALADARES TAVARES 
MÁQUINA DE TURING PARA CONVERTER EXPRESSÕES 
EM NOTAÇÃO «INFIX» EM EXPRESSÕES EM NOTAÇÃO 


«POSTFIX» — II. ATRIBUIÇÃO DE VARIÁVEIS DE ESTADO 
ENTRADA E SAÍDA. 


Turing Machine for the translation from «infix» to «postfix» 
notation — II. State, input and output assignments. 
L. T. MAGALHÃES 
M. P. O. RICOU 
O TRANSISTOR -—- REACTOR ATÓMICO., 
The transistor — The atomic reactor. 
BENJAMIN DA CONCEIÇÃO MENDONÇA 
"THE PERFORMANCE OF ACOUSTO-OPTIC CELLS FOR LA- 
SER-DOPPLER ANEMOMETRY. 


Performance de células acusto-ópticas para a anemometria com 
Laser. 

D. F. G. DURÃO 
Jd. H. WHITELAW 


Edições da Técnica 


TABELA DE PREÇOS 


TABELAS PARA CÁLCULO DE BETÃO ARMADO 


J, S. Brazão Farinha, 6.º edição — 1970 , 375800 
TOPOGRAFIA GERAL 
C. Xerez— 1.º vol, 2.º edição— 1970 . . (*] 
C. Xerez— 2.º vol. 2.º edição — 1966 . 180800 
TABELAS TÉCNICAS 
Vicente e Eng.º Brasão Farinha, 8.º edição — 

1977 E") 
MANUAL DE HIDRÁULICA GERAL 
A. Lencastre, 2.º edição — 1969 375800 
DICIONÁRIO DE UNIDADES E TABELAS 
DE CONVERSÃO 
Vasco Costa e Osvaldo Francês, 1.º edição — 

1959 35$00 


(*) LIVROS A PUBLICAR BREVEMENTE 


Desconto de 10% aos assinantes 


GUIA DE ANALISE QUÍMICA DAS ÁGUAS 


A. Herculano de Carvalho, 1.º edição — 1961 90$00 
CALCULO DE PÓRTICOS 
G. Kani — Trad. da 7.º edição alemã — 1962 75$00 
MANUAL DE ESTRUTURAS 
Brasão Farinha — 1.º VOLUME 
Encadernado 625$00 
Plastificado h87$50 
Ze O dd VOMNÕE mz uúrrgws a (*) 
BETÃO PREESFORÇADO 
(A cargo dos participantes na Semana de Be- 

tão Preesforçado — |. S. T.— 1975) ... (*) 


PEDIDOS À [EENIRO 


NUMERO 438 FEVEREIRO DE 1977 ANO Lil VOLUME XXXVIH PUBLICAÇÃO MENSAL 


DIRECTOR SUMÁRIO 


Delmar Baptista 


COLABORADORES 
283 — M. ABREU FARO — Estabelecimento e interpre- 


António Maria da Fonseca tação física da transformação de Lorentz. 


José Agostinho Marques Physical establishment and interpretation of the Lorentz 


Jorga Braz transformation. 
Armando Ruano 


Manuel M. S. Anastácio 
José Luis S. V. Azevedo 
António Ricardo Pinto Gonçalves 


295 — ANTÓNIO GOUVÊA PORTELA — Introdução 
aos conjuntos vagos. 


Introduction to fuzzy sets. 
Henrique José Carvalhão 


Sérgio Duarte Santos (F. C. Tec- 
nologia —- Coimbra) 


311 — ERNST BECKER — Education and research in 
mechanics in the Federal Republic of Germany. 


Ensino e investigação em mecânica na República Fede- 


DIRECÇÃO REDACÇÃO 
ral Alemã. 


E ADMINISTRAÇÃO 


Av. Rovisco Pais, |. S. T. — Lisboa 


Telefone 72 93 23 317 — LUÍS VALADARES TAVARES — Continuous hy- 


drological time-eries discretization. 


PROPRIETÁRIO | Discretização de séries hidrológicas continuas. 
Ma Es du Bi Ti | 
| Ê 
| 325 — L. T. MAGALHÃES e M. P. O. RICOU — Má- 
x | quina de Turing para converter expressões em 
ASSINATURAS 5 nº 10 nº | notação «infix» em expressões em notação 
be h. n. 


«postfix» — Il. Atribuição de variáveis de esta- 
do, entrada e saida. 

Turing Machine for the translation from «infix» to «post- 
fix» notation — Il. State, input and output assignments. 


Continente e Ilhas 200$00 3508. 
Países de lingua 
Portuguesa e Es- 


panha 2208600 38058. 
Estrangeiro ...... ao 4505. 

| 335 — BENJAMIM DA CONCEIÇÃO MENDONÇA — 
Número avulso...  — 608. O transistor — Reactor atómico. 


Não sa publica em Agosto E The transistor — The atomic reactor. 


Setembro, 
359 — D. F. G. DURÃO e J. H. WHITELAW — The 


performance of acousto-optic cells for laser- 
-doppler anemometry. 


Performance de células acusto-Ópticas para a anemo- 
metria com Laser. 


Os artigos assinados são da 
exclusiva responsabilidade dos au- 
tores. 


de 


COMPOSIÇÃO E IMPRESSÃO 


COOP. ESPIRITO SANTO, SCARL 
Rua do Espirito Santo, 19-A - Odivelas 


a 


INDICE 


Fanafel 
MAGUE 


Novobra, Ida. . 


ALFABÉTICO 


Pág 


IV 


Il 


IV 


DOS ANUNCIANTES 


Sociedade Portuguesa Cavan . 
SM. PP; LOA. , nes ú 
Sondagens e Fundações A. Cavaco 
sopecate . 


Pág. 


II 
II 


US 


A SKF não tem apenas rolamentos de esferas 


e de rolos qa o) 1) ) p 
Temos também rolamentos de agulhas E y 

rótulas Ty e embutes: 5). 

Evidentemente, temas do mesmo modo toda a espécie 


OM E 


Além disso a SKF duna, muitos outros produtos 


de acessórios 


para rolamentos. 


de qualidade. 


Rodas e carretos conicos, por exemplo & 


Fusos roscados de rolamento Wifccaiii 


precisao que transformam o movimento de rotação em 


linear. 


Fabricamos também machos 
e brocas sESG= . 
O mesmo se passa com ponto ivos =>. 
e passa c pontos rotativos gy gsm 
E muitos sabem que também temos um bom Serviço 


Técnico. 


SOCIEDADE SS LIMITADA 


LISBOA = PRAÇA DA ALEGRIA, 66-A PORTO - RUA DELFIM FERREIRA, 604 
TELEF.: 362301 - TELEGR.: ESKAEF - TELEX:|2/56 TELEF.: 692054 - TELEGR.: ESKAEF 


TÉCNICA | 


Mi ÁGUAS 
k SUBTERRÂNEAS 


ESTUDOS - PESQUISAS 
CAPTAÇÕES 


para colaborar 
Di na reconstrução 


A o Pais 


EMA te 


" | o 
Ebal PRA E 
A r da mi = e 

, “im a Z 
"A MI |. nn 


? j pe 
| ' mo: Wo. 
SONDAGENS E FUNDAÇÕES A.CAVACO,LDA. 


O e 
RUA RODRIGO DA FONSECA, 62 : LISBOA - 1 E 


TEL. 56 11 71/4 - TELEX 12 476 ACAVAC P* TELEG. ACAVACO 


A ic E fr Ma 
| // 
4///) my S 


Postes de cimento 


armado Caran 


Qualidade que 
desafia o tempo 


Av. Visconde Valinor: 76-1.º - Tel. 766014 (7 linhas) Lisboa 
TÉCNICA II 


* FUNDAÇÕES ESPECIAIS 
- ESTUDOS GEOTÉCNICOS 
= CAPTAÇÕES DE ÁGUA 


Sopecate 


ns ESPECIALISTA DESDE 1947 


“elit R.DO ARSENAL, 146 -2' 
cia Telf. 320208-364010-360437 


Wed” LISBOA-2 


Já pensou. 


«-. Que O seu anúncio 


é muito eficaz se for 
publicado na revista que 
os seus clientes lâem? 

Os milhares da exemplares 


de cada número da nossa revista | 


são lidos por engenheiros 

8 técnicos, por pessoas ligadas 
à indústria, à investigação, ao 
ensino, pelos seus clientes, 

em suma. Consulte-nos 


ENGENHARIA 


Avenida Rovisco Pais - Lisboa 1 - Telef, 129323 | 


NUMERO 438 


FEVEREIRO DE 1977 


ANO Lil 


VOLUME XXXVI 


PUBLICAÇÃO MENSAL 


C.D.U. 530.12:513.75/77 


Estabelecimento e Interpretação Física da 


Transformação de Lorentz 


SUMÁRIO 


Procede-se neste trabalho ao estabelecimento e 
interpretação fisica da Transformação de Lorentz. Para 
tanto, recorremos aos dois postulados fundamentais da 
Relatividade Restrita e aos conceitos de velocidade e 
dessincronização. 

Na base do método utilizado, a Transformação de 
Lorentz ganha significação fisica o que permite uma 
análise directa do seu formalismo e, decorrentemente, 
a sua escrita imediata a partir do significado alcançado. 

Este é o segundo de três artigos que estamos a 
publicar na Técnica sobre a Transformação de Lorentz, 
sua interpretação fisica, e métodos geométricos para a 
determinação da referida transformação. 


1 — INTRODUÇÃO 


A Transformação de Lorentz pode e tem sido esta- 
belecida a partir de diversos métodos, 


A cada método corresponde, em geral, uma deter- 
minada atitude, e, muitas vezes, uma finalidade bem 
definida. 


Tratando-se da transformação de coordenadas que 
presido à Relatividade Restrita, é natural que de entre 
esses métodos existam aqueles que se filiam directa- 
menta nos dois postulados da Relatividade Restrita: 


— Postulado de um Princípio de Relatividade Res- 
trita, que traduz o facto experimental de as leis da Física 
serem as mesmas, em conteúdo e forma, para diversos 
referenciais mas restrito, na sua extensão, aos referen- 
ciais de inércia: qualquer um onde se verifique a Lei 
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SUMARY 


In this paper the establishment and the physical 
interpretation of the Lorentz Transformation is underta- 
ken. To this end we have but used concepts of velocity 
and desynchronization and the two basic postulates of 
Special Relativity. 

On the basis of the method used, the Lorentz 
Transformation gains physical significance, which allows 
a straightforward analysis of its formalism and, conse- 
quently, its immediate writing following the meaning 
arrived at. 

This is the second of the three articles we are 
publishing in Técnica on the Lorentz Transformation, its 
physical interpretation and geometric methods to deter- 
mine ft. 


da Inércia e todos aqueles que em relação ao primeiro 
estejam em movimento rectilineo e uniforme, 


— Postulado da Invariância da Velocidade de Pro- 
pagação da Luz no Vácuo, que traduz o facto experi- 
mental de a velocidade de propagação de qualquer sinal 
electromagnético, no vácuo, ser uma constante universal 
“c” para todos os referenciais de inércia. 


Utilizaremos como fonte primeira os dois postula- 
dos da Relatividade Restrita. 

Na Relatividade o espaço tem carácter relativo e 
apenas faz sentido falar do espaço relativo a um dado 
referencial S onde, um dado acontecimento, fica carac- 
terizado por quatro coordenadas E(X,Y,Z,T). Na 
Relatividade Restrita a geometria euclidiana é a que 
conduz à maior simplicidade de trato e tem a vantagem 
de nos manter conceitos da geometria dos corpos sóli- 
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dos. Por estes factos se adoptam coordenadas cartesi- 
anas ortogonais, melhor diriamos, euclidianas se em vez 
de T se considerar cT. 


Reconhece-se que não é possível por uma experiên- 
cia local individualizar qualquer ponto ou direcção do 
aspaço relativo de um referencial pelo que o continuum 
aspaço-tempo é aceite como homogéneo e isotrópico, 

Posto isto, consideremos dois referenciais de inér- 
cia stx,yz.t) e S(X,Y,Z,T) em que o primeiro, 
a na perspectiva do segundo, se acha animado de mo- 
vimento rectilíneo e uniforme de velocidade v. 

Admitamos que um dado acontecimento E é assl- 
nalado no espaço relativo do referencial 5 pelas coor- 
denadas 


Es 1X YZ T) 


Essa acontecimento também pertencerá ao espaço 
relativo de s onde é assinado por 


E lx y z t) 
O problema fundamental consiste na determinação 
da lei de transformação 
X= A (x y. z t) 
Y Y My 21) 


£ [XY € 1) (1) 


II 


Z 


I 


T=Tlixy z t) 


que permita para quaisquer valores (x, y, z, t) das coor- 
denadas de um acontecimento, em s, determinar as 
coordenadas (X, Y, Z, T) do mesmo acontecimento, 
em 5. 

E também, como é óbvio, resolver o problema 
inverso. 


Pelo primeiro postulado a transformação (1) deve 
ser linear, 


Por outro lado e uma vez que o espaço-tempo é 
homogéneo não há perda de generalidada em arbitrar 
quaisquer pontos de S ou s como origens das coor- 
denadas desses referenciais. 

Sendo assim, podemos impor que as origens das 
coordenadas dos dois referenciais, s e S, assinalem 
um mesmo conhecimento. 

Isto significa que, na transformação em causa, a 


RÃ=Y=Z 


| 
Õ 
D 
= 
| 
Õ 


corresponde 
Xx=y=z=o0 e 


Arbitrando assim, a transformação, que já se vlu 
ser linear, é também homogénea. 

Uma vez que o espaço relativo de qualquer refe- 
rencial é isotrópico, a orientação dos eixos pode ser 
qualquer e assim admitiremos que para T=t=o os 
eixos espaciais dos dois referenciais estão apostos e 
dirigidos segundo os mesmos sentidos, respectivamente. 

Em seguida, podemos orientar o eixo dos XX de 
modo a que a origem das coordenadas espaciais de s, 
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x=y=z=o, seja vista mover-se sobre o eixo dos XX 
vom a velocidade v. 

A partir daqui, recorrendo uma vez mais ao primeiro 
postulado e isotropia do espaço, é fácil demonstrar que 
a transformação pretendida se reduz à forma: 


X = Ax + Bt 

T = Bx + CG 

Y=Y (2) 
Z=azI 


A lei de transformação que resulta para as coor- 
denadas reveste uma forma especial directamente ligada 
à orientação relativa que se adoptou para os eixos e 
origem das coordenadas arbitradas, por isso se diz 
Transformação Especial de Lorentz. 

será a partir de (2) que actuaremos e os nossos 
objectivos fundamentais são dois: 


— Estabelecimento da Transformação de Lorentz 


— Interpretação física da Transformação de Lorentz, 
a partir do método utilizado para o seu estabelecimento 
e dos conceitos de velocidade relativa e dessincroniza- 
ção relativa, 

Atingidos os objectivos visados, apresentaremos 
um exemplo onde se esclarece quanto anteriormente 
sa expôs. 


2 — ESTABELECIMENTO DA TRANSFORMAÇÃO ES- 
PECIAL DE LORENTZ 


Como se disse, com generalidada e por simples 
recurso aos dois postulados, é possível determinar com- 
pletamente a Transformação de Lorentz. 

Para atingir os objectivos que visamos utilizaremos 
um procedimento diverso. 

Assim, abandonamos a generalidade para passar- 
mos a considerar uma situação concreta da qual resul- 
tará uma transformação univoca de coordenadas. Uma 
vez que nessa situação se actuará em conformidade com 
os dois postulados, a transformação obtida é, como seria 
da esperar, a Transformação Especial de Lorentz, uma 
vez que o estabelecimento é válido para quaisquer valo- 
res das coordenadas, 

Antes de mais e como resulta de (2) a transforma- 
ção especial deixa invariantes as coordenadas ortogonais 
à direcção da velocidade v e resume-se à forma simples: 


X = Ax + Bt (3 a) 


T = Ex + Dt (3 b) 


que correlaciona entre si as coordenadas segundo a 
direcção de v e os tempos te T dos dois referenciais. 

Vê-se, desde já, que o tempo perde o carácter 
absoluto que tinha na Transformação de Galileu. 

Recorrendo a uma onda electromagnética plana po- 
demos, num dado referencial, definir geometricamente 
acontecimentos simultâneos. Isto é, traçar lugares geo- 
métricos de acontecimentos simultâneos e de tempo T 
em relação à origem das coordenadas: 
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Amv=ZeoaTE os, 


Adivinhando-se dificuldade no estabelecimento 
dessa onda e sua utilização, poderíamos emitir no ins- 
tante T=o, e a partir de um plano espacial de refe- 
rência, sinais electromagnéticos segundo a normal posi- 
tiva do plano. 

Esses sinais existem num plano que, com este 
esclarecimento e a partir daqui, designaremos por frente 
de onda. 

Consideremos a fig.- 1 onde, em relação ao refe- 
rencial S (X, Y, Z, T), se ilustra a definição geométrica 
de lugares de T = const. 

Na figura indicam-se, de entre outros, dois proce- 
dimentos possíveis. No primeiro, utiliza-se uma onda 
plana com a normal de onda dirigida segundo o sentido 
positivo do eixo do YY. Então, e de acordo com o 
princípio de invariância da velocidade de propagação 
dos sinais electromagnéticos, a frente de onda, lugar 
geométrico dos sinais emitidos, caminha a partir de 
Y =o e durante o intervalo de tempo 


FIG. 1 — Dois procesoss de estabelecer lugares geométricos de 
tempo constante a partir de sinais electromagnéticos 


E 
E D 

uma distância 
ct, 


Para T, = o essa superfície de onda coincidia 
com o plano de referência Y = 0 
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Sendo assim, um acontecimento E, localizado no 
plano XY, Z = o, tem por coordenadas 


E Nati 
(4) 


E,. por exemplo, é um acontecimento de coordenadas 
Ea (Xo 0, 0, 0) 


uma vez que corresponde a T, = O 

Para o que se pretende basta-nos o plano XY pelo 
que passaremos a considerar a intersecção da superfícia 
de onda com o plano Z = o. Isto significa que a recta 


Y = cT, 
Z=o0 


traduz um lugar geométrico de acontecimentos ocorrendo 
no tempo T, e localizados no plano XY. 

No segundo procedimento utiliza-se uma onda pla- 
na com a normal dirigida segundo a direcção a. 

Nesta representação, E, é um acontecimento de 
coordenadas 


E, (X,. Ejs O, Tt; 


cT, — cos « X, (5) 
Y = > cm mes 
” sen « 


Y, obtém-se facilmente a partir da projecção de OX,E, 
sobre o eixo dos TT. 

As cordenadas do acontecimento E, correspondente 
a T;=o serão, de acordo com (5), 


E, (Ho Yo 050) 


dm Xo cos a sa Xo 
sen a tan e 

Fica deste modo estabelecido um processo de 
obter lugares geométricos de tempo constante, a partir 
do recurso a «ondas planasa. 

Uma vez que a partir daqui apenas usaremos o 
plano dos XY no referencial S e o plano dos xy no 
referencial s as frentes de onda serão linhas rectas, 

Posto isto, estamos em condições de actuar em 
relação à seguinte situação: 

No referencial s(x,y, z, t) estabelece-se uma 
onda plana com a normal de onda dirigida segundo 
a direcção e sentido positivo do eixo dos yy. Mais 
concretamente, uma linha recta que para t= O coin- 
cide com y=0, 


Para t=0, a frente de onda que se estabeleceu 
acende, na origem das coordenadas desse referencial, 
uma lâmpada que permanece acesa. 


Ao fim de um tempo t essa frente de onda atingiu 
o plano y=ct e na sua progressão vai acendendo 
lâmpadas colocadas uniformemente ao longo do eixo 
dos yy, lâmpadas que após uma breve emissão de luz 
se apagam. 
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A lâmpada colocada na origem das coordenadas 
cor esponde uma sucessão continua de acontecimentos 
de coordenadas 


E, [o,0,0,t] 


e em que o acender da lâmpada corresponda a t=o0. 
Esta lâmpada é pois uma marca vwespacial» da 
origem das coordenadas, x= y=z=o, que resta as- 
sinalada para todo o tempo t>o, 
Ao acender de cada uma das outras lâmpadas, 
colocadas no eixo dos yy correspondem acontecimentos 
de coordenadas 


E, [o,ct,0,1] 


A cada t corresponde uma lâmpada, aquela que se 
situa à distância y = ct da origem. 

Teremos, neste caso, um sinal luminoso que per- 
corre o eixo dos yy com velocidade relativa c, velo- 
cidade da frente de onda que comanda o acender de 
cada uma das lâmpadas. 

Estabeleceu-se assim no eixo dos yy um indicador 
do tempo t do referencial s, tempo relativo à origem 
dos tempos t=o. 

Nos termos do primeiro processo descrito e repre- 
sentado na fig. 1, as rectas y = ct são lugares geomé- 
tricos de acontecimentos de tempo constante. 

De acordo com o segundo postulado, a onda elec- 
tromagnética estabelecida, em s, é também, em 5, uma 
onda electromagnética propagando-se com a velocidade 
c relativa a este referencial, 

Na fig. 2, no alto e à esquerda, representa-se o 
que em S se observará: 

De acordo com o conceito de velocidade relativa 
a lâmpada que assinala a origem «espacial» de s deslo- 
ca-se ao longo de X com velocidade v. 

De acordo com o segundo postulado o sinal lumi- 
noso correspondente às lâmpadas situadas em y, é ob- 
servado em S como um sinal luminoso com velocidade 
«c» que durante o tempo T se deslocou a partir da 
origem de uma distância cT a segundo uma direcção 
a, à determinar, 

Comparemos as coordenadas dos dois aconteci- 
mentos. De acordo com o que se disse: 

A lâmpada na origem x=y=zZ=0, é no tempo 
t, corresponde E. (0,0,0,t). 

À lâmpada em x=272=0, em y=ect, e no tempo 
t, corresponde E. [o, ct, o, t]. 

Chega-se assim à conclusão que estes aconteci- 
mentos têm iguais coordenadas, (x) = 0 e (t) = t, pelo 
que, de acordo com a forma particular da Transforma- 
ção Especial: 


X = Ax + Bt (3 a) 


T = Cx+ Dt (3 b) 


lhes correspondem, em S, acontecimentos com Iguais 
KH eoeT. 

Para o tempo T, o em S, à lâmpada que no ins- 
tante t se encontra acesa em Y = ct, corresponde um 
acontecimento afastado cT da origem, segundo uma 
direcção «, e com uma abcissa X = vT, uma vez que 
esta é abcissa da lâmpada na origem, x=y =2=0,6 
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se demonstrou que estes dois acontecimentos têm 
iguais coordenadas X e T. 

É o que se representa, no alto e à esquerda na 
fig. 2, e daí se conclui que 


V 


Cos au = 
c 

A partir daqui podemos utilizar a onda electro- 
magnética estabelecida em s, e tal como é vista em 5, 
para definir lugares geométricos de tempo constante T 
nos termos do segundo procedimento representado na 
fig. 1. 

Esses lugares geométricos serão rectas ortogonais 
à direcção «, o que resulta da conjugação do segundo 
postulado com o facto das transformações serem li- 
neares, 

Por outro lado, e uma vez que na transformação 
em causa as distâncias segundo direcções ortogonais 
à da velocidade ficam invariantes, segue-se que uma 
recta xy, z= o, lugar geométrico de tempo constante 
y=ct é observada em S segundo uma recta em XY, 
Z=o, de equação 


Y=y=et 
E | w E" [O 
- * * 
É, 5 
=) =] 
» * 


xa, 
FIG. 2 — Estabelecimento da Transformação de Lorentz a partir 


do conceito de velocidade relativa e dos dois postulados da 
Relatividade Restrita 


TÉCNICA 438 


Mas sendo assim, só haverá um valor para a 
coordenada X de dois acontecimentos com tempos t 
e T nos referenciais s e S, respectivamente, 

Tudo isto se representa no segundo desenho, a 
contar do alto e à esquerda, da fig 2. 

A determinação de X é imediata, basta projectar 
a linha quebrada OXE sobre o eixo dos TT, obtendo-se: 


cT=X cosa+ct sena 
de onde resulta: 
cT — ct sena 


cos a (6) 
em que 


v 


COS a =— 


Na fig. 2, ao alto e à direita, representa-se outra 
experiência, idêntica à que se descreveu. 

Nesta, a onda electromagnética é estabelecida em 
S com a no:mal de onda dirigida segundo o eixo dos 
YY aproveitada depois para definir lugares geomé- 
tricos de tempo constante em s. 

Seguindo a mesma via chegava-se à conclusão 
que a variável x ficava determinada pela relação 


ct-—cT sena 
-= COS e 
em que (7) 


COS a = 


Esta relação seria de esperar e está em directa 
concordância com o Princípio de Relatividade. 

Do facto, na primeira experiência era s que, na 
perspectiva de S, estava animado de movimento rela- 
tivo e uniforme de velocidadea v. 

Pelo primeiro postulado é igualmente legitimo 
dizer-se que, na perspectiva de s, é S que está ani- 
mado de movimento relativo e uniforme de velocidade 
[— v). 

Sendo assim, seria de esperar que a uma relação 


X=FMT,t, v) 
em S, correspondesse, em s, uma relação 
x=F lt, T, —v) 
Note-se que na segunda experiência o ângulo entre as 


direcções e sentidos positivos de x e té a=T-a 
pelo que: 


vz 
sen « = sen rt 7 — 
Cc? 


cos e == cosa=-—- v/c 


É o cosseno, função impar de v, que troca de sinal 
quando v passa a (— v). 
A partir das relações obtidas: 
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cT—- ct sena 


x = 
cos « (8 a) 
ct— cl sena 
agia — cos e (8 b) 


imediatamente se obtém: 


x + ct cosa 


+ 
sen « (9 a) 
ct+x cosa 
5 a E 
sen a (9 b) 
como se obteria: 
X— cT cosa 
di e 
sen a (10 a) 
cT —- X cose 
CE = —em em 
sen a (10 b) 


Exprimindo sena e cosa em v/c 


V 


Cos a = 


RE Sá (11) 
sena = 1 ma a — 
cº 


poderemos ainda escrever: 


A = coca ds (x + vt) ax 
Vime 2 a 
1 V 
T= (t + a x) 
V 1 — v2 /c2 (12 b) 
1 
X= aa - (X— vT) 
V 1 — v?/c? (13 a) 
/ 
1 V 
t= SN —-—— 4] 
C- (13 b) 


| / 4- v2/c? 


que, com as relações, 
Y=y 
Z=z 


traduzem, como é bem conhecido, a Transformação 
Especial de Lorentz. 
Obtém-se, deste modo, o resultado que se esperava. 
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Prosseguindo, e em directa relação com a fig. 2, 
verifica-se que, rodando a representação correspondente 
à segunda experiência, as duas representações se ca- 
sam consistentemente e conduzem à figura, em baixo 
e ao centro, onde o acontecimento 


E 


de tempos relativos t e T abcissa relativa X tem uma 
representação, na primeira experiência, que coincide 
com a representação associada à segunda experiência 
em que o mesmo acontecimento de tempos relativos t 
e T é assinalado na abcissa relativa x. 

Da reunião das duas figuras resulta a sobreposição 
dos lugares geométricos, da t e T constantes, associa- 
dos a cada uma das experiências, pelo que fica esta- 
belecido um processo geométrico de obter as trans- 
formações 


[x, t] & [X, 7] 


Na fig. 3, mostram-se as representações que resul- 
tam para a resolução dos seguintes problemas 


FIG. 3 —- Geometria associada ao método utilizado. Ressalta de 
qualquer das figuras, a Invariância da Forma Quadrática 
Fundamental 


— Conhecidos [X, T] determinar [x, t] 
— Conhecidos [x, t] determinar [X, T] 


É traduz a direcção a associada à 1.º experiência 
(m — «) traduz a direcção (m — a) associada à 2.º ex- 
periência, 

Destas figuras resulta imediatamente: 


— Da primeira: 
oc Tt+x=Xº + ce? t? 
— Da segunda: 
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ettx=x + Tº 
São relações idênticas que associadas a 
+ My 
Z=2z 


traduzem a invariância da Foima Quadrática Funda- 
mental: 


(xe +y+z)—-c te=(Xº+Y2+ Zº) —e T? (15) 


Posto isto estamos em condições de passar à inter- 
pretação física da Transformação de Lorentz. 


3 — INTERPRETAÇÃO FÍSICA DA TRANSFOR- 
MAÇÃO DE LORENTZ 


Consideremos a fig. 2 onde se representam as 
duas experiências que nos conduziram à Transformação 
Especial de Lorentz. 

Como se salientou recorreu-se, então, ao conceito 
de velocidade relativa o que, em última análise, sig- 
nífica dizer que no referencial S a origem das coorde- 
nadas espaciais de s se desloca com velocidade v, me- 
dida através de adequados registos em X e em T. 

Este conceito é, em certa medida, o que assiste 
à Transformação de Galileu. Está certo mas não é tudo. 

Na realidade, pode pôr-se questão semelhante em 
relação à origem da coordenada tempo do referencial s. 

Para analisar este assunto, admitamos que no 
referencial S uma recta de tempo constante comanda 
sobre o eixo dos TT o sobra o eixo dos XX o acender 
de lâmpadas que lá se imaginam dispostas e com den- 
sidades tais que podemos admitir que no eixo dos TT 
e correspondente ponto do eixo dos XX existem lâm- 
padas que se acendem simultaneamente sempre que a 
frenta de onda por elas passa. 

Atendendo a que 


cosa= v/c 


imediatamente se conclui da figura que essa frente de 
onda, ver a fig. 2, põe em correspondência a abcissa X 
com um tempo 

v 


EE— 


cº 


Deste modo, e de acordo com a definição de rectas 
de tempo, o primei:o sinal luminoso, sucessão de lâm- 
padas que se acendem, progride no eixo dos TT com 
velocidade c. 

Em contrapartida o sinal luminoso no eixo dos XX 
progride com velocidade 


X/ 


Esta velocidade, maior do que c, provém de lâm- 
padas independentes que se acendem sob o comando 
de uma frente de onda que se desloca com velocidade 
c segundo a sua normal de onda e que pela inclinação 
u« deste com o eixo dos XX corta este eixo em pontos 
que se deslocam com a velocidade 


TÉCNICA 438 


Cc E Fi 
cos a v V 
C 


Trata-se de um artifício em relação a fontes inde- 
pendentes que não contraria, portanto, o segundo pos- 
tulado. 

Posto isto, notemos que o eixo dos XX coincide 
com uma recta de tempo t constante. Mais precisa- 
mente, t= 0. 

Isto significa que ao acende: da lâmpada na ab- 
cissa X corresponde, em s, um acontecimento de co- 
ordenadas 

E |x,0,0,0) 


Em S esse mssmo acontecimento tem por coor- 


denadas 
v 
Es( X,0,0, ee) 


Isto é, para x= 0 et=o tinhamos 
Xk=0 e To 
enquanto que para x = 0 a t= o temos 


V 


X&Eo aoT= X 


Cc? 
Conclui-se daqui que a origem do tempo de s, 
t=o, se separa de T=o,e ao longo de X. Essa sepa- 


ração, medida em S, é dada por: 


V 


T-0sT= X 


cê 


Definiremos como dessincronização a variação de 
tempo por unidade de comprimento. 

A dessincronização relativa da origem do tempo 
de s, t= o, será então 


v v 
d="—XX= 
cf Cc (16) 


e tem, como é evidente, as dimensões do inverso de 
uma velocidade. 

Para a outra experiência, S em movimento e s 
fixo, seria o mesmo, apenas se observaria uma mudança 
de sinal para a dessicronização que passaria a ser 


V 


c? (17) 


Do que se disse podemos inferir, em última aná- 
lise, o seguinte: 


— É tão legitimo dizer que um referencial está 
animado de velocidade relativa v em relação a outro, 
como dizer que está animado de dessincronização rela- 


V 
tiva d= E em relação a esse outro referencial. 
EC 
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— Na Transformação de Galileu e para a veloci- 
dade relativa v não havia a contrapartida da dessin- 
cronização relativa d. Na realidade, a dessincronização 
era nula o que conferia carácter absoluto ao tempo. 


—— Na Transformação de Lorentz os dois conceitos 
coexistem a só por experiências particulares que fixem 
um ponto do espaço x = x, ou um instante do tempo 
t=t, se pode evidenciar, respectivamente, uma velo- 
cidade relativa pura e uma dessincronização relativa 
pura. 


— Da um modo geral a origem das coordenadas de 
5 separa-se, no espaço e no tempo, da origem das co- 
ordenadas de S: 

No espaço, com velocidade v ao longo da tempo T. 

No tempo, com dessincronização d ao longo do 
espaço X. 


— Sendo assim, é de esperar que a um aconteci- 
mento E com coordenadas, em s. 


E, (x, t) 
corresponda um acontecimento, em 5, com coo:denadas 
Es (X, T) 


em que X seja dado por 


R=SvTÃHK,. x (18 a) 
a T dado por 
W 
TEA Ret (18 b) 


—. vT traduz o que em X, e ao longo de T, a origem 
x= o se deslocou. 


w ' 
e M traduz o que, em T, e ao longo da X, 
E 
a origem t = o se separou no tempo. 
— K,x e Kt representam a contribuição respecti- 
vamente para X e T das coordenadas x e t. 


—— Na Transformação de Galileu K. = K, =1 
e v/cº = o; tudo se passava como se existisse uma via 
fisica de conciliar os dois referenciais e a que corres- 
pondesse c=&, 


—— Na Transformação de Lorentz e dada a simili- 
tude de trato que confere ao espaço e ao tempo será 


de esperar que 
vê 
c? ) (19) 


De —— , dada a isotropia do espaço-tempo. 


Lil 


—— Tudo isto está claramente evidenciado na fig. 4 
de onde se conclui, imediatamente, que: 
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FIG. 4 — Interpretação da Transformação de Lorentz a partir dos 
conceitos de velocidade relativa e «dessincronização relativa 


VARE: 
K=kK, =K, =sena = ia Ea (20) 


— Este último resultado encerra, em si, e por 
conveniente interpretação, os conhecidos efeitos relati- 
vísticos da «contracção das dimensões segundo a di- 
recção do movimento» e edilatação do tempo», mas 
disso não nos ocuparemos aqui. 


—— Apenas nos interessa salientar que às separa- 
ções no espaço e no tempo há que juntar as projecções 
sob:e os respectivo eixo, e na orientação implícita no 
método que se utilizou, da coordenada x e da coorde- 
nada t. 


— Por essa projecção as coordenadas x e t refle- 
tem-se em X e T através do mesmo coeficiente 


Conclusão: 
Tendo em consideração o exposto e com a signi- 


ficação adquirida, as Transformações de Lorentz escre- 
vem-se imediatamente: 


x=v+ ave x 
v 

Pre MA) | =Neig q 
c? / 


8 reciprocamente 


(21 a) 


[21 b) 
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x=-—vt+ VÁ — v2fc? X 


(22 a) 
v Di caca 
[= ——— 4 N q-vi/c? T 
c? | (22 b) 
De (22 a) e (22 b) resulta: 
1 
Xx= ———— (x + vt) 
a — ville? 
Are Ra (23 a) 
1 v 
E CR 
/ = (tt) 
= (23 b) 
De (21 a) e (21 b) resulta 
1 
= ——— (X — vT] 
Vi-S 
1 (24 6) 
1 V 
t=>—————— (T- — x 
Wide ni 
| — e (24 b) 


que são as formas usuais que a Transformação Especial 
da Lorentz reveste. 


4 — EXEMPLO 


Consideramos de interesse ilustrar com um exem- 
plo aquilo que se expôs relativamente ao estabeleci- 
mento da Transformação de Lorentz e sua interpretação 
física. 

No alto da fig. 5 representa-se, no referencial 
s, as seguintes sequências de acontecimentos: 


— Nos pontos 1, 2, 3 localizadas sobre o eixo 
dos xx acenderam-se, durante um intervalo de tempo 
muito curto, três lâmpadas. 

Acendem-so simultaneamente no instante t=o. 


— Seguidamente, a lâmpada situada em 3 torna 
a acende:-se no instante t, e uma vez mais no instante 
=2t,. 


— Para qualquer destes acontecimentos y=o 
pois as lâmpadas localizam-se no eixo dos xx. Como 
as representações de y = o coincidem com as det=o 
esses pontos têm um quadrado sobreposto. 


-—— Uma vez que nos situamos no plano dos xy 
será, para qualquer acontecimento, Z = o. 


-—— Portanto, as coordenadas dos acontecimentos 
referidos são: 
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E, [0,0,0,0] E, [X,,0,0,0] E, [X,. 0, 0,0] 


Es [Xs, O, O, t,1 Es [xs O, D, t,] 


— Além destes, outros acontecimentos se passa- 
ram. Assim, uma frente de onda de normal dirigida 
segundo o eixo dos yy vai progredindo com velocidade 
c e acendendo nos instantes t,, t,, t, etc., lâmpadas 
situadas no eixo dos yy. 

As unidades de espaço e tempo são tais que 
conduzem a c= 1, pelo que os pontos figurativos de 
y coincidem com os de t. Esse facto vai assinalando, 
também, por um quadrado, 

— Sando assim, rectas ortogonais ao eixo dos yy 
traduzem lugares geométricos de tempo constante: 


Po: 3 passa a recta y=o 
Por 4 passa a recta y = ct, 


Por 5 passa a recta y = ct, 


Posto isto consideremos o referencial S onde se 
observa o referencial s em movimento rectilíneo e uni- 
forme de velocidade v = 0,6 c. 

Aplicando a Transformação Especial de Lorentz 
obtêm-se facilmetne as coordenadas desses aconteci- 
mentos no referencial S. É o que se representa na parte 
inferior da fig. 5. 

O desenho vai feito à escala e é fácil comprovar 
a situação dos pontos, ou seja as suas coordenadas. 

Analisemos agora o que se observa em 5: 


— Ao longo de uma linha quebrada e designada 
por T indicam-se as coordenadas (X, T) dos aconteci- 
mentos E,, E,, E, E, e E, em 8. 


— Ao longo do eixo dos X e numa sequência de 
estrelas brancas assinalam-se as coordenadas (X, Y) 
dos mesmos acontecimentos. 


— A sequência T, T, T, traduz as coordenadas 
(X, T) dos acontecimentos ocorrendo em Y,. Yar Yys 


— À sequência Y, Y, Y, mostra-nos os locais 
onde as lâmpadas colocadas em y, y, Y, Se vêem acen- 
der em S, trata-se de uma representação em (X, Y). 


— bs rectas ortogonais ao eixo dos T, ou Y tra- 
duzem rectas de tempo constante T. As coordenadas 
dessas rectas, que não se representaram, são cT,, CT, 
CT,- 


— Mas de acordo com o segundo procedimento 
considerado na fig-1, a recta dos Y, Y, Y, pode cons- 
tituir-se num eixo dos tempos T, direcção AB, e as 
rectas AD e BE, ortogonais a AB, são então, e tam- 
bém, rectas da tempo constante T., 


— Note-se que as distâncias: 
OY, OY; OYq o. 


coincidem com as ordenadas 


DÃs, Ole CE 
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— Note-se que OBEDO é exactamente a relação 
geométrica que aparece na fig. 3, ao alto, e mais por- 
meno:rizadamente na fig. 4. 


FIG. 5 — Sequências de acontecimentos e sua representação 
coordenada nos referenciais s e 5 


— Sendo assim, e em relação ao acontecimen- 
to E,: 

OB = cT, ordenada de 5 no plano (X,T) 

OD= X 

ED = ct 


BE ==x abcissa de 3, em s, e também de 4 ou 5 


V 
OA = X 
Cc: 
OC= nvT 


—— Assim se verifica claramente a teoria do mé- 
todo que sa apresentou, 


Finalmente, interpretamos à luz do que se expôs, a 
fenomenologia observada em 5. 
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— Três acontecimentos em s, t= o, são vistos 
em S, nos tempos assinalados pelas ordenadas de 1, 
2 e 3, mais sugestivamente assinalados nos relógios. 


— A velocidade do sinal luminoso segundo X 


é 1/d 
cs 


V 
É o fenómeno puro da dessincronização relativa, 


—. Três acontecimentos em s, e com x= o, têm 
as suas abcissas em S deslocando-se com velocidade 
v. abcissas de Y, Y, Ya. 

É o fenómeno puro da velocidade relativa, 


— À sequência 3, 4, 5, em s, e que traduz variação 
no tempo, em s, para x = conste. corresponde uma 
dilatação do tempo, em S, através de um factor 


17 V17— v2/c2. 


— A distância BE = x, projecta-se em S e na sua 
contribuição para x, através de um factor de redução 


o 


— O tempo DE = ct, projecta-se em S e na sua 
contribuição para T, através de um factor de redução 


va 
j — — 


= K, mt : 
Ge 


— Estes dois factos foram devidamente assinala- 
dos na interpretação da Transformação de Lorentz. 

A figura, devidamente interpretada, revela ainda 
outros factos mas julgamos que estes que se referiram 
são os mais importantes. 


5 — CONCLUSÃO 


— A partir da interpretação feita tem significado 
escrever directamente: 


=== 
=vT EE 1— v2/c? X 


V 


T= X + V = vio ? 


c2 

— Interpretando assim, é fácil escrever a Transfor- 
mação de Lorentz para outra orientação de eixos ou 
outras origens das coordenadas. Bastará atender ao 
significado intrínseco da transformação em que X e x 
representam distâncias medidas segundo a direcção da 
velocidade, v, e T e t rep:esentam os tempos dos refe- 
renciais 5 e s. 


— às distâncias ortogonais à direcção da veloci- 
dade, v, ficam invariantes. 


292 


—— Na referida escrita está implícita a separação 
vT que a Transformação de Galileu contemplava, e 


v 
figura, e adquice naturalidade, a separação a 


—— À estas separações estão associados dois con- 
ceitos: o de velocidade e o de dessincronização. 


— Para v= o obtemos a transformação unidade 
NX =x 
T=t 


o que, pela velocidade ser nula, traduz a imobilidade 
relativa no espaço e, de acordo com a geometria dos 
sólidos, que 


x =x 


Também, e pelo facto da dessincronização ser nula, 
traduz que 


Te à 


o que significa um tempo comum aos dois referenciais 
e encerra, em si, o conceito de simultaneidade, só 
aplicável a acontecimentos do mesmo referencial, con- 
ceito a partir do qual faz sentido falar do tempo de 
um referencial e da existência de relógios síncronos que, 
em cada lugar, o indicam. 


-—— No âmbito de uma Teoria de Sistemas 


v 


E ;— 
V ja 

cr 
traduz, em certa medida, um factor físico da ligação. 


Quando v=c os dois sistemas ficam desligados 
e, para qualquer x ou t, teremos: 


R=oeT 
X 
T= — 
E 


o que significa que x e t têm projecção nula em 5. 
A relação 


x = et 


condiciona as coordenadas da origem de s, em 5. 

Assim, os acontecimentos do espaço relativo de s 
passam todos a situar-se, em 5, em planos ortogonais 
ao eixo dos XX, satisfazendo as suas coordenadas às 
relações: 


X=cT, para qualquer x ou t 
X 

T= —— , para qualquer x ou t 
C 

Y=y 

Z=z 
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Evidentemente, que outras interpretações são pos- 
síveis, mas estas, que nos parecem fundamentais, são 
aquelas que nos vêm da significação que emerge da 
Transformação de Lorentz quando escrita: 


A=vT + Va -ure X 
ú e 
q + Vavo é 
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"Introdução aos conjuntos vagos”? (continuação) 


CAPÍTULO 5 
RELAÇÕES BINÁRIAS 


5.1 DEFINIÇÃO E PROPRIEDADES 

Se X-Y = então XX Y =iWxAibt=L? e É? será 
o conjunto referencial. De novo (2, v,4A) seja um 
reticulado e J a família das aplicações funcionais de ll? 
em il. 


As principais propriedades que podem possuir as 
Relações Binárias são: 


5.1.1 SIMETRIA 


So (u,.U,) ES 62 Ja (u, u,) = Ja (u,,u,) então 
A é uma relação simétrica. 


5.1.2 REFLEXIVIDADE 

Sel (uu) & W for J, (u,u) = 1 então A é uma 
relação reflexiva (repare-se que não basta ser 
Ja (u, L) > O). 


5.1.3 TRANSITIVIDADE (DIRECTA) 


Se Ú(u,.u,), (u,.u,) E WU? for 


v | Ja UU, Us) | MEO (u,. u, ) dp fu; Ê Us) | 
| 
= Ja (4; Us) 


então A é uma relação transitiva. 
Note-se que a transitividade se testa pela compo- 
sição de À o A; com efeito 


AoAã= Vv 1go!! ss A 
o af am u; Jo Ja tu; vo |] 


5.1.4 TRANSITIVIDADE (INVERSA) 


Sel (u,. U,), (u, u)E WU? for: 


António Gouvêa Portela 
Prof. Catedrático |. S. T. 


| Ja lu, us). a ROO ) Jalu; e) | 


= Ja (U,, Us) 
então À é transitiva (inversa). 
Esta definição só pode ser aplicada se o conjunto 


reticulado for complementado ou, quando muito, pseudo- 
-complementado (à Zadeh). 


Os resultados das duas formas, directa e inversa, 
não são, em geral, os mesmos. 


5.1.5  ANTI-SIMETRIA 


So (uu) & for: 
E U,) 55 Jp (Us vu | U. 


U E tu, U,) = Ja lu, U,) | =0 


então A é anti-simétrica. 


5.1.6  ANTI-SIMETRIA A ZADEH 


Só se aplica a reticulados ((0, 1), *,v). 


Sel (u,, u,) E 6 for 


- 


| Ja (u,, 5) > o | = Ê (Us, U,) = o| 
então À é anti-simétrica à Zadeh. 
5.1.7 FECHO TRANSITIVO DUMA RELAÇÃO BINÁRIA 


Seja A C qQf e simbolize-se por ALA=A 0 A. O 


elemento genérico Ja, (uj.u;)S A? terá a forma: 


(*) Continuação do trabalho «Introdução aos Conjuntos Vagos» publicado na Técnica-433, Fevereiro de 1976, pág. 201 a 217. 
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Jas (u,, U;) E AC = 


= v 


u, | s| JalupeUy do Jal e) | 
= = 


e se À for transitiva diecta já vimos que Ja V Jas = 
= Jyouseja J,; > Jjy ou À2 Aº. (veja-sa 4.1.3). 


Dum modo geral A* 5 ATI. 


Chama-se fecho transitivo duma relação binária a 


ni a 
à = 


a) Se A for transitivo então A é transitivo. 
t 
b) Se A 2A então À é transitivo. 
t 


c) Se A = A“ então À é transitivo. 


d) Se a partir da ordem k, Art = ak então 


k 
A= UA 


são propriedades que resultam da definição imediata- 
mente, 


5.1.8 CAMINHOS NUM GRAFO FINITO 


Seja dado: 


Um grafo AC qb”, cujo caracterizante Ja, CE (OQ, v,4) 
reticulado, onde Card dé é finito. 

Define-se caminho finito Cj de u, para u, , num 
grafo finito A, um aq-uplo, com a seguinte forma e 
propriedades: 


a) cu =(u,,... uu, ) é um q-uplo ordenado, 
cujos elementos vão ser resimbolizados, passando a 
escrever-se: 


E — (u, =. 


q u) 


com ju=u, e qu=U, (g >2 e finito). 


q 
bVje(t, que 


Cc) Us. pe(i ... q), ju o pu 
donde resultará que 


Sd lijriem.q=1): 


onde O é o Ínfimo de 2. 
O comprimento h de c, é dado (q — 1): 
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A largura | do caminho c,. é definida da forma 


seguinte: 


| (Cj) E [Jal ps qui] comjell..g— 1) Ie 


9.1.9 CAMINHO DE MAIOR LARGURA DE u, PARA u,.. 


Seja E ne o conjunto vulgar de todos os caminhos 


distintos iniciados em u e terminados em Uj. + 


: Tm SÊ ne tf : 7 
Os caminhos Ci. E= O o Cj = É a dizem-sea 


distintos se os respectivos q-uplos forem distintos, isto é, 
tiverem pelo menos um elemento distinto, ou a ordem 
é diferente, ou q é diferente. Então o caminho de 


maior largura será o caminho Cl tal que: 


e ttcjo. tejo | = teto 


L. 


e A | | lei). ci) | = ei) 


Ve 8 “e A AS 
para Y Cie E og O GT Sik 


Simboliza-se | (cj) = Li. 


5.1.10 CAMINHOS DE MAIOR LARGURA E DE UM 
DADO COMPRIMENTO h. 


O conjunto Fido todos os caminhos iniciados em 
a 


u, e acabados em u, pode ser dividido em classes 
de equivalência usando o comp:rimento h como critério 
de classificação, 


Assim Du representará todos os caminhos inicia- 


dos em U, e terminados em up que tenham um com- 
primento h e 


2 = 
nO = Oi 


Pode definir-se a largura , Li. 
maior largura e de comprimento h 


h 
do caminho da 


nbik = ci) 


que satisfaz ao seguinte conjunto de condições: 


(O 
a) Ci E h É a 
b) v | ci), tejo) | Hei) 
e A | Hei), ci) | = ci) 


CO 
4 pa — PA PE 1 
para Ácii E nO CF Ci 
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5.1.11 APLICAÇÃO DESTES CONCEITOS A GRAFOS 
CUJOS RETICULADOS (£2, vw, 4) ASSOCIADOS SÃO 
TRANSITIVOS, COMO SUCEDE COM OS RETICULA- 
DOS DE ZADEH.* 


A transitividade exprime-se pela proposição: 


Va, pjes? | (sta 8) = a). (ato, B)= 8) |es «<p 


Então teremos: 


a) Ile) = 4 | Sad ue qu) | [Etta eg 1 
, - 


Haverá pelo menos um Ja(,u;py U) = cj) 
com LE, a g—1) 


b) Lig =Hejg) = v Hei) Mei) )= Nei) 


=: á 
v Cik “= te 


é suficiente para definir L,, 
c) Há pelo menos um JA ( mL mijU=L,, 
ME sa g =) 


d) porque 


Le =Y I A | Ja GussqU) ]| 


e da definição de fecho transitivo (5.4.3) pode estabele- 


cer-se que: 
Jahlujuç) =, Li (directamente) 
e) Ja (up up )=La (directamente) 

f) Se h> Card qf, então 


para m « Card Qf. Basta lembrar que o grafo é finito 
e que caminhos com h >n têem repetições (circula- 
ções) que não alteram o resultado devido à idempotên- 
cia dos operadores 4 e V., 


n 
g) is () A en=cCard q 
i=1 


pelas razões já invocadas, 


(*) Veja-se Anexo 2. 
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5.2 — DEFINIÇÕES DE VÁRIOS TIPOS DE RELAÇÕES 


5.2.1 PRE-ORDEM VAGA 


Pre-ordem vaga é uma relação binária, transitiva e 
reflexiva. Daqui resultam as seguintes propriedades: 


a) Aº=AoàA e 


— 


Ja” (uu) = 4 | Sa (sou, DA dp (o | 


1 À. 
mas para i = 1 será 
Ja lu u)áSo (u,u,) 
e porque Ja (u,,u,) = 1 f(reflexividade) será 
Ja lu u,) A Ja lu u,)=J, (u,u,) 


Atendendo à transitividade 


Ja lu, U,) > = 
donde 
Ja? (uu) =Ja (u,u,) 
e finalmente 


A: = A (o não apenas Aº* c- A) 
b) Dum modo geral AE =A e A=A (feho tran- 


sitivo) t 


5.2.2 SEMI PRE-ORDEM VAGA 


Relação binária não reflexiva e transitiva. 


5.2.3 PRÉ-ORDEM ANTI-REFLEXIVA VAGA 


Transitiva mas fue 6: Ja lu,u)=0 


5.2.4 RELAÇÃO DE SIMILITUDE (OU SEMELHANÇA) 
VAGA 


É uma relação p:é-ordem vaga simétrica. 
Porque é simétrica pode escrever-se 
«=Ja (lu u)=Ja (uu) 
B=Ja tusug)=Ja (us, u,) 
vY=Ja luu)=Ja (ug u,) 
Dadas as relações de transitividade será: 
a>zB AY 
B>ae Ay 
Y>2«4 B 


donde o seguinte conjunto de restrições: 
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a > B B>]eNrNvY>2cD «= 
azBAB2cN TABS A=B 
a TABDTA TADRSA=B= 
C>7nN827Y A TAB V= 


a 

> 

2YNB>De A TV2AD7E=A 
> 
- 


YABar À 
Y Ne Sã AA 


a Y20e50=7 
> 


PB > B=7 


pi 


Estas propriedades resultam da simetria e podem resu- 
mir-se; 


(la>B=W)VIBDe=y)ViyDa= ep) 


5.2.5 SUB-RELAÇÕES DE SIMILITUDE DUMA PRÉ- 
“ORDEM VAGA 


Seja dada uma p:é-ordem vaga À € q(º e é, c Ab 
um sub-conjunto vulgar que constitui um sub-referencial 
e J,€C J um sub-conjunto de J e relativado a Qt; se 


Tx tuu=3 a lu) Vuu,e U, 


os elementos de “(., formarão uma sub-relação de simi- 


lidade. 
Uma sub-relação de similituda é maximal se não 
estiver contida noutra sub-relação de similitude. 


A classe de todas as sub-relações de similitude 
maximais se existir numa pré-ordem designa-se por 
classe de similitudes da prá-ordem. 


Uma pré-ordem que pode ser decomponível em 
classes de similitude diz-se uma pré-ordem redutível. 
5.2.6  ANTI-SIMETRIA (VAGA) 

Viu,.u) CA qu <u: 
[Ja fu,u,)=5 plus, U,)] ou [J fu, u,) = Ja (us, u,) = 0] 
diz-se que a relação à anti-simétrica, 


5.2.7 GRAFO ANTI-SIMÉTRICO VULGAR ASSOCIA- 


DO A UMA RELAÇÃO VAGA 
Vu, u,) E WU? 
THU, E Us | Ja (up us) > Ja (ug u,) | =» 
[ui u,) & Gelu, u;) é G | 
2) u, = Us | Ja (u,, Ug) = Ja (us, uy) | = 


> | (u,.U,) e (u,,U,) é G | 
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5.2.8 ANTI-SIMETRIA PERFEITA (VAGA) 


Define-se Vu, Us) E Té? 8 (u, + U,) 


será: se J, (u,,u,) >0 entãoJ, (u,,u,) = O 


5.2.9 RELAÇÃO DE ORDEM VAGA 


É uma relação binária reflexiva, transitiva, anti- 
simétrica. Porque J, (u,.U,) >J q (u,s.u,) em Q% in- 
duz-se uma ordem: 


U >U; 


5.2.10 RELAÇÃO DE ORDEM TOTAL VAGA 


Se o grafo vulgar associado representa uma or- 
dem total. 


5.2.11 RELAÇÃO DE ORDEM PARCIAL VAGA 


Basta que o grafo vulgar associado seja de ordem 
parcial, 


5.2.12 RELAÇÃO DE ORDEM PERFEITA [VAGA) 


Se se usar para definição de anti-simetria a de- 
finição 5.2.8. 


5.2.13 RELAÇÃO DE ORDEM ESTRICTA (VAGA) 


transitiva 
anti-reflexiva 


Relação estricta | 
anti-simétrica 


transitiva 
reflexiva 
anti-simétrica 


Relação não estricta 


5.3 — OUTRA RELAÇÕES TÍPICAS 


5.3.1 RELAÇÕES ORDINAIS (VAGAS) 


São relações que gozam das seguintes proprie- 
dades; 

— reflexiva 

— anti-simétrica (veja-se 5.2-6, 7, 8) 


— não tem circuitos no grafo vulgar associado de 
comprimento superior a 1. 


5.3.2 FUNÇÃO ORDINAL DUM GRAFO VAGO 


Veja-se Anexo-3 que trata este tema para grafos 
vulgares. 


A definição de relação de ordem vaga ordinal im- 
põe que goze das seguintes propriedades: 


—— reflexiva 
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— anti-simétrica 

— sem circuitos 

— transitiva 

É possível definir uma função ordinal para a rela- 
ção de ordem vaga ordinal. Para o efeito usa-se o 


método descrito no Anexo-3 para onde se remete O 
leitor interessado nesse método. 


5.3.3 RELAÇÃO DE DISSEMELHANÇA (OU DISSIMI- 
LITUDE) VAGA 


Caracteriza-se por ser: 
— simétrica 
— gnti-transitiva 


Falta definir o que se entende por anti-transitiva. 
Anti-transitividade define-se como segue: 


| (u, V), (v, ww), (u, W') = Ai? 


Ja tum) < 4 [Jalev) vd) iv w)] 
v 


5.3.4. Tem interesse em confrontar a «dissemelhança» 
com a «semelhança» 


q =Jaluw)=Ja (w, uU) 


Semelhança: simetria B=.alu, v) = Ja tv, u) 
e Ja lv w) -— Ja (we, V) 

| o a= Ja (lu w)=Jy (wu) 
Dissemelhança simetria ( g = Ja tu, v) = Já (v, u) 


v=Jalvw)=J; (wv) 


São ambas simétricas. 


Dissemelhança Jalu,u)=V VAQ=-O 
Semelhança Jã (uu) =O0 Vvy O=V 
Finalmente 
Semelhança ARS 
ses O a a transitividade 
Dissemelhanço «S 4 (B, VIy |) 
v 


a aplicação das leis de De Morgan verificam as ex- 
pressões de transitividade. 

Portanto se o reticulado for complementado e se 
Ja satisfizer as regras de semelhança então Jq ve- 
rifica as regras de dissemelhança, 

Pode estender-se esta afirmação aos reticulados 
de Zadeh ([0,1], 7.%) que gozam da propriedade da 
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pseudo complementariedade. À demonstração seguia os 
mesmos passos. 


5.3.5 Relação parecença é reflexiva a simétrica. 


5.3.6 Relação de não parecença é anti-relexiva e si- 
métrica. 


5.4 — RELAÇÕES E SUAS PROPRIEDADES QUANDO 
O RETICULADO É DE ZADEH ([0,1], 1,1%) 


Embora se tenham faito algumas referências ao 
reticulado de Zadeh e porque este é dos mais empre- 
gados em conjuntos v2gos, convém dedicar uma alinea 
ao estudo de algumas situações típicas que ocorrem 
com as relações já apresentadas no caso do reticulado 
ser de Zadeh. 


5.4.1. — CONCEITO DE DISTÂNCIA MIN MAX (D) 
ENTRE ELEMENTOS DUMA RELAÇÃO BINÁRIA. 


daluv)=1—Ja (uv) 
com uv eU e Ja E ([0:1), £ ,%) 
e se À fo: o complemento de A 

da luv) =Jz (u,v) 

Estas definições serão aplicadas a vários tipos de 
relações e há que provar, em cada caso, se a definição 
de da satisfaz as regras de uma distância. 

Aplicação a uma relação de semelhança (ver) 

dalu v) = 1 — Jalu,v) 


sendo À uma relação de semelhança vaga 
A é: transitiva (directa) 


reflexiva 


simétrica 
Porque é transitiva teremos: 
Jatuv) > y DD, tu w)4Jalw, v)] 
ww 


ou 


1— Jalu,v) <4 [1 Ju w)yp (1 Ja (w,v))] 
w 
ou 


da (uv) < 4 [da lu w) pda tw, v)] 
ww 
Porque é simétrica: 


Ja lu, v) = Jalv,u) 
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donde 


daluv)=1—Jatuv)=1—Ja(vu)=dalvu) 


Porque é reflexiva: 


Ja tu,u) = 1 (supremo) 
donde 
datuu) =1—J tu =1—1=0 (ínfimo) 
e da=0 5 1—Jaluw)=0 ou Jaluw)=1 


e daí que uUu=w 


logo é mesmo uma distância. 
Finalmente tem interesse em reparar que se desig- 


narmos por A a relação complementa; de À com J, = 
=1— Ji então Ja lu, v) = Jx (u,v) e ainda que o 
operador que conjuga d, (u, w) com da (w, v) 6: 


A, [da lu, we) by dalw, u)] 
e daí a designação «min. máx». 


Aplicação a uma relação de parecença (ver 5.3.5). 

Uma relação de parecença A é reflexiva e simé- 
trica; daí que o fecho transitivo A, é ainda transitiva 
(directa) portanto A, é simétrica reflexiva e daí uma 
relação de similitude (ou semelhança). 

Define-se distância entre dois elementos u,v de 
A, como sendo 


da luv) =1 — dida 


5.4.2 — CONCEITO DE DISTÂNCIA MINSOMA D EN- 
TRE ELEMENTOS DUMA RELAÇÃO BINÁRIA. 


Os operadores são agora “ e = (soma disjuntiva] 
Se A for o conjunto dado seja A* o conjunto ob- 
tido de A pela cperação seguinte: 


Sae luv) =%7 [a tu w) E Joe lw v)] 
wW 


que designaremos por min soma. 


Aplicação a uma relação de parecença 


D 4 ju, v]j=1 — né (u, vw) 
Ja, (u, v) — “A, Iv, u) > DA (u, v) 
Ja, (uyu)=15 Da, (uu) =O0 


a finalmenta 


Ja luv) > y Ja, (u,w). Ja (w,v)] 
t t 


Wi 
Ou 
T=5% tuv)>r =, luw) (1, (mw v)] 
t ui t t 
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1-7 tuvM> vit (Ia lu w)+Js (wv)- 
t Nú t t 


— a, (u, Wi] dj, tw, v))] 


Sã <+t DG tuw) EJA (wv)] 
t w t t 


como se desejava demonstrar 
D (uv) >D(uw) 6) D (w,v) 
onde =) é a ope.ação minsoma no domínio do reticulado 


de Zadah. 


5.4.3 — FECHOS TRANSITIVOS 


a) Fecho transitivo do tipo máx. prod. 


Os operadores definidos no reticulado são (y,*º) 
ou seja Máx. e prod. vulgares. 

Se A for um conjunto vago, numa relação a com- 
posição do tipo máx, prod. será: 


Jas — JaoA = y [JA (u, w). Ja (we, v)] 
Wy 


o se Card A=y então o fecho transitivo (max. prod.) 
será: 


Outras formas de compor típicas são: 


Ja?= J(A0A) =p [Ja tu v)éJda tw v)] 
ng 
ou 
Ja = J(AOA) — 17 [JA tu, w) 4 Ja (Wo, v)] 
wW 
Note-se que: 


Ja (vw). Ja (wmv) < Ja tu, wW) Ja (av, V) 
e daí que, se for verdade Aº c A pela operação (W.% ) 
então também pela operação (y,*) será verdade que 
AI A. 

Também se pode concluir que 

(At >2A)D(A<O >A) 

onde <y,* > e <yº> são os operado:es max-min 
e max-prod. respectivamente, 


5.44 — COMPOSIÇÃO E DECOMPOSIÇÃO 


Aplicação destas operações a certos tipos da rela- 
ções binárias. 


a) Relação de similitude 
Dum modo geral será: 


=7va.R, 
[ra 
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ondeoSasS1l ea.R, é o producto vulgar a pelos 
elementos de R, .R. é o conjunto vulgar de nível a. 

Demonstrar que se R é uma relação de similitude 
será; 


(a, o &s) rea R x, — (Ra,) 


e que, dum modo geral, R, é uma relação de equiva- 
lência. Para demonstrar que R, é uma relação de equi- 
valência basta verificar se é reflexiva simétrica e tran- 
sitiva o que é fácil sabendo que uma relação de simi- 
litude é também reflexiva simétrica e transitiva, 6 que 


a, > Cc Ra DRa, 
é igualmente imediata a demonstração. 
b) Relação de ordem perfeita. 
É fácil ver que se 


q > % > Ra, — Ra, 


c) Aplicação das composições ao caso em que 
4 = [a;b]JcR. 
Seja 


Jalxy)=a aefl 


Jply=p pen 


sdejem y e 4 as cperações definidas em £! e o operador 
de composição, por exemplo, A 


à ed) ESA 


J(a,B) (x, z) — [JS à (x, y) à Jp ly, z)] 
y 


e 
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CA ant am 
Sejam P, P, P, os pontos de cruzamento; 
há que estudar [0,P,] [P,. P.] [P,.V] 
Em cada intervalo tomamos: 
edi O E da it 


ada ti) > sda MM 
as partições [P,;Ps ] são função de (x, z) 


Ja (9) sa 


ou ,Jp (lv) se 


ANEXO 1 
A1.1 — Correspondências (relações) 

Sejam dados dois ronjuntos JL, e Al, e forme-se 
o producto cartesiano Q(, X l,; seja GC QL, X ks, 
um sub-conjunto dado; então G é o grafo da correspon- 
dência IT estabelecida do conjunto ql, para o con- 
junto Q... 


A1.2 — Aplicação é toda a correspondência que gosa 
da propriedade adicional de: 


Vxeuwlyvel,:lxy)EGe GC, X UW, 


A1.3 — Uma aplicação diz-se funcional ou função se 
a só se existir apenas um y para cada x, no grafo G: 


(xew ye qu: (xy) EG 


A1.4 — No texto, para simplicidade de linguagem o 
vocábulo aplicação vai ser usado por aplicação funcional. 


A1.5 — Relação binária, diz-se das relações cujos gra- 
fos G estão contidos no producto cartesiano IL XxX É 
que simbolizaremos por É. 


A1.6 — Conjuntos Ordenados. 


Uma relação binária w definida em 6? diz-se re- 
lação de ordem se for transitiva e equivalente, 


a) Diz-se transitiva se: 


(xy) e (y,z) E w então (x,z) E w 
b) Diz-se equivalente se: 
(x,y) e (y,z) E w então x=y 
A1.7 — Pré-ordem 
Uma relação binária diz-se uma pré-ordem se for 


reflexiva o transitiva, 
A reflexividade define-se como: 


(xeU xx) E q 
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A1.8 — A noção de inclusão ( € ) permite ordenar o 


conjunto das partes de É, ED (AE). 
Com efeito, sendo dados A e B E |D(u) como 
á possível definir BC A ou À E B será então possivel 


ordenar os elementos de 9 (SE). 
A1.9 — Relações de ordem mais usadas e seu signi- 
ficado semântico, 


a) (>) —x é maior ou iguala y:x =, y 


b) (>)-x é maior que y:x >y 
Assim x >y=(x>Dy)Vlix=y) 


4110 — Cofinal (coinicial) 


Num conjunto ordenado (qt, >) ao conjunto 
X c qi dá-se o nome de cofinal (ou coinicial) se se 
verificar: 


Vyeqixex:ix>y (ouy >s) 


A1,11 — Elemento minimal (maximal) em (UM, >) 


O elemento a & q/,é minimal (ou maximal) se, 
sendo a > x (ou x > a), então tivermos x = a 


A1.12 — Elemento menor (maior) 


ae (9é >>) será elemento menor (ou maior) 


is, 


sol xe tor x > a (ou x< a). 


Ã1.13 — Conjunto minorante [majorante ) 


O elemento x(q(., >) é minorante (ou majoran- 


te) do conjunto À E dl se y EAfory >x (0ux>y). 
O conjunto minorante m (ou majorante M) de A 
é o conjunto contituído por todos os elementos x tais 


que Y YE À tem-se y >x (ou x > y). 


A1.14 — Limite inferior (superior) 
Limite inferior (superior) é o elemento maio; (me- 


nor) do conjunto minoranta (majorante). 


41.15 — Conjunto filtrante 


Todo o conjunto A, pré-ordenado, é filtrante à 
direita (à esquerda) se toda e qualquer parte finita de 
A não vazia for majorada (minorada). 


(*) Outras designações para reticulado Vemos, por exemplo: 


— LÁATIS (ou LATA) expressão portuguesa 
— REDE ORDENADA, 
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A1.16 — Reticulado (*) é qualquer conjunto ordenado 
Ti. onde todo o À C 9X tem limite superior e limite 
inferior, e esses limites pertencem a “k. 


ANEXO 2 — CONJUNTOS RETICULADOS 


A 2,1 — Definição 


Seja dada uma triada (!), v, 4) onde Q é um con- 
junto com dois ou mais elementos e v e À são duas 
cparações internas conjugadas definidas em 1? a que 
gozam das seguintes propriedades: 


a) fechadas e unicidade no resultado 
Ula p)er 
Ylapjen 


vila B)=- 5E80 
Ale, B)=y€0 


5 e y são únicos. 


b) vv, A são operações associativas 
vLlevtB,y)] =v Lv (a B),7] 
Aafe MB, =ATA (a B),7] 


para Va p,yel. 

c) vw, à são operações comutativas 
Vila B)=viB,a) 
A (a, 8) = A (B,o) 
para (a. Ber 


d) Existência e unicidade dum supremo e dum 
infimo. 

Existo um só elemento VE l tal que 

Vaco,vlaV)=V (supremo) 
doen:heg,a(4,0)=0 (íntimo) 
e) lIdempotentes 


V (mae)=a 


Ale a)=a 


para | qEen 


ft) Absorventes 
A [av la Bi]=ea 


7 [Ale Bl] 
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Para Ta, Pe t 


Em geral o conceito semântico destes simbolos 
operatórios são respectivamente: 


4 — limite inferior 


w — limite superior 


A 2.2 — Propriedades mais interessantes dos reticulados 


a) Em consequência da comutatividade e da asso- 
ciatividado é fácil ver o seguinte: 


AtaTa(a,8),7],6) =p pode escrever-se sim- 
bolicamente 4 (x, 2,7, 8) = p tendo esta representação 
um significado claro, 

Se ainda designarmos por A = [e 8,7, 9) e AC 
então 

ATA] =», 


ese B=[aw,y,m.0) com B € 92 então podemos escre- 
var expressões tais como: 
ATiA NBJU (y)] 


O mesmo se diz para o operador q. Em qualquer 
dos casos há que determinar os elementos correspon- 
dentes e que satisfazem à expressão entre os pa;êntesis 
rectos antes de operar com v ou á, 


b) Sejam 
ABES 


8,71,9 EM 


v [B]=0 
v EA =% 
v [AUB]= 8 
Então é fácil ver que 
v (9,76) =8 
podendo eventualmente 9 e (ou) y serem iguais a 8. O 
mesmo cem à, 
c) Supremo e infimo 
Seja: 
a EN AACS, CardÃ>2 
Então teremos: 
vIAV)=V 
StA,V) = A (A) 
v (AO) =v (A) 


a (A,0)=0 
v tla,0)=« 
Ata, 0)=0 
v le V)=V 
A(eV)=a 
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d) Relações de ordem e pre-ordem parciais 


Vamos estabelecer a seguinte equivalência entre 
as operações internas (vw), (4) e a relação de or- 
dem (<). 


sendo (e, B) E £º” seja 


[Ale B)=a)UlyiaB)=8)] C5> las B] 


Está assim definido um nexo entre os referidos 
símbolos operatórios e a relação de ordem. 

Só quando for verdade o primeiro membro é pos- 
sível definir uma relação de ordem e nem para todos 
Os pares (a, B)E 2,4 (a, B) e/ou v (a, B) são « e/ou 
B: em muitos casos v l(e,B)=5 e/ou Ale, B)=Y 
sendo 8 Ea e BB y%a e B. Então entre (a, B) não se 
pode estabelecer uma relação de ordem: daí que esta 
ordem seja parcial. 

Mais adiante demonstrar-se-à a transitividade, equi- 
valência e reflexividade da relação (<) e ainda que 
ela não é total em geral, mas o seu grafo não é vazio. 

Antes de passar a essas demonstrações é fácil 
de ve:ificar a verdade das seguintes expressões: 


[A (e B)=8] 5 
[vleB)=7] 5 


[(5 <a) n(ô < 8)] 
[a <Yyn(B<7] 
Com efeito: 

ô=AleB)=A(A (ee), B)= 


=A(AleB)a)=A(5,a) 
mas 


[4 (85,0)=5]€5 [5 < e] 
Repetia-se igual raciocínio para 8 e obtinha-se [5 <. 8] 
g.e. d, 
O mesmo para a operação (v ). 


Vejamos agora se o operado: < é efectivamente 
uma relação de ordem parcial e uma pré-ordem parcial. 


Com efeito 
|—. Tiansitividade 
Se for («<B) e (B<y) então (a < 7) 
Com efeito 
la< B)€> US (a, B)=a)U ly (a, 8) = 8)] 
(A<NESUA (By) =B)U ly (8,7) =7)] 
Há as quatro seguintes hipóteses a estudar: 
[4 (a, 8) =] N TA (8,7) = 8] 
[4 (a, 8) = 0] NALv(B.7) = 7] 
[vla B)= BIN TA(B.y) = 8] 
[v le 8) =8]NTv (8.7) = 7] 
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Vejamos para a primeira: 


Ate y)=4 [Ala B),y] = 4 [0,4 14B,7)] = 
=áAlefBj=e 


[4 (a, 7) = a] 
Para as formas v (8,7) =%y basta recordar que 


ATB,vIBy)l=B8 


ou 
AB y)=8 
[viBev)=7) 5 AMB.7)=B 


Pode acontecer que nem ála, 8)= a nemy (a, 8) = 
=B e então não se pode declarar que (a < A). 
Isto significa que haverá pares (a, 8) que não perten- 
cem ao grafo da relação de ordem e daí que a ordem 
não é total, em geral, 

Por outro lado existirão sempre um V e um O 
tais que 


VlaV)=V 
(ren 
à te, 0) = O 
e pelo menos se pode-á declarar que 
u<NV 
| ae! 
O<e 


daí a existência de uma ordem parcial, 


|| — Equivalência 


Esta propriadada expressa-se pela seguinte frase 


[la<B)N(B<a]ES [a =] 
Com efeito 
[Ate B)=)UlviaeB) =8)] ESla< 8) 
[A (e, B)=B)U (via) =] ES (Pp <e) 
le<BIN(p<a) E> 
€> [la (e B)=0)N (A (a, B)=B)U 
Ullate,B)=acn (via B)=a)]U 
Ulivia,B)=B)nlAleB)=8)]U 
Ulivia B)=BIntiviaB)= al] 


| 


ou 
(la=8)Ula=8)Ula=8)U(a=ãçB) 
Logo 
(e=8) > [les B)n (BS al] 


HI — Reflexividade 
Vlaaj=e 


VeceL,Alma=e a 


Dal que (<) seja também uma pré-ordem. A pré- 
-ordem e a ordem definidas são parciais, como já se 
disse. 


304 


e) Conceito de máximo e minimo 

Quando for: 

Vila gjem falegp)=e 
. R (a, B) = B 


via B)=B 
V te B)=a 


então os símbolos A e vw serão substituídos por, res- 
pectivamente, + e y e correspondem ao conceito se- 
mântico «máximo» e «minimo», 

Note-se que se 8 e at la, B) e jogando com as 
propriedades da associatividade e comutatividade, te- 
mos que: 


A(Aj)=5 580 (em geral 5€A) 
mas 
A(A]=:e sE ACO 
Identicamente 
vYiA)j=7 yeAhfEO 


com Card A >2 

Pode ainda escrever-se: 
"AM =0 
v O) =» 


ai =0 
em vez de | 
v(O)=V 


A.2.3 — Propriedades adicionais dos operadores A e q. 


Ã.2.3.14 — Distributividade 
Os reticulados dizem-se distributivos se: 
v [lr A(B.y)]=á4[vtapB)v (e 7) 
e Ale viBy)]=v lálaB),4 (a, 7)] 


quaisquer que sejam «,/f,Yy pertencentes a 2, 
O grafo seguinte não é distributivo 


V 


2.3.2 — Complementariedade 
B é complementar de a se for 
V le B)=V 
e S (ae, B)=0 
Se f for único para cada a então £ diz-se o 


complementar de a e representa-se por a. 

Caso esta propriedade e respectiva unicidade seja 
verdadeira para todos os elementos de 2, (2, w,4) 
diz-se complementado. 
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A principal propriedade daqui resultante é a da 
involução: 


(a)=a V 20EU 
Com efeito 
v la, e) =V 


A le,e)=0 (por definição) 


donde 

v La), 0] =V 

A La), (a)] = O 
mas o complementar de a é único, logo 


(a) = a Vae 2 


A2.3.3. — Pseudo-complementariedade 


Sea «sn for um elemento relacionado com «EL 


de tal modo que (o) =e, isto é, a aplicação sucessiva 
do operador (—) reproduz o elemento de partida (ou 
seja, (—) goza da propriedade da involução), então 
diz-se que se pode definir em (2, v, A) uma pseudo- 
-complementaridade. 

É evidente que se não exige que 


v lma)=V e A la 0) =0 


A2.3.4. — Teoremas de Morgan 


Os teoremas apresentam-se da forma seguinte: 
Teorema 1— v (a, B)=A (a, B) 
Teorema 2— A je, B) = Va B) 


Para que se verifiquem os teoremas de De Morgan 
é necessário que: 


— O reticulado seja complementado. 


Va E 92, exista um só elemento B =ae em 
tal que: 
v lae)=V 
Far (e, a) =0 
— Se designarmos por 
vila Bl=y, yE28 
AleB)=5, 8€EQR 


como (£24,w, à) é complementado y e ô existem e são 
únicos. 


—— Finalmente 
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p é y são igualmente únicos. 
Em resumo: 


V le, p)jef? existirão e serão únicos: e, B, Y.8, 
p, w ao todo oito elementos contando com a e f. Estes 
oito elementos estão entre si vinculados por oito pro- 
posições: 


v lwa)=V 
A (a,a) = 0 
vi B)=V 
A(B.B)=0 
vivia B), ale B)l=V o vIiypl=V 
alvo), ala B)]J=O0 o A[yp]=0 
vIAteB), vlapl=V o vIS&vyvI=V 
AlAteB), via B)]J=0 o A[8,4]=0 


Não há muitos reticulados que satisfaçam às con- 
dições de Da Morgan. Como exemplos teremos: 


V V 
V 
0 0 0 


fig. 1 fig. 2 fig. 3 


Por simples verificação é fácil de ver que são 
complementados a que as leis de De Morgan se aplicam. 

O reticulado da figura 1 é o reticulado de Bool. 

Os teoremas de De Morgan também se aplicam 
a certos reticulados pseudo-complementados, como o 
reticulado de Zadeh assim definido: 


([0,1],W, é ), distributivo e pseudo-complementado. 


Este reticulado goza das propriedades de lhe ser 
aplicável os teoremas de De Morgan. 
Com efeito, definindo 


Vaen:ig=1-a E] (a) = « 


(=) símbolo do pseudo-complemento. 


então: 


| 


Vl Bj=y(1-041-B)= 


II 


= à (a, B) 


o mesmo se podendo fazer em relação a (a, 6). 
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A2.3.5 — Reticulado Modular 


Se Va, BYE Q for: 
(la<B)> aviBgáy)=levB)Ay 


diz-se que (º2, 7, 4) é um reticulado modular. 
Como exemplo temos: 


A2.3.6 — Conceito de dominância e comparabilidade 


n 
Seja dado o produto cartesiano > A, e sejam 
[= 
(jo Xge cer Xp), 


igualmente dados dois elementos 


(jo Vigo cuco Yá ) onde x; 8 y; EIA, Y e (1... n). 


Diz-se que (x,,....x, ) domina (y,,....Yy ) se: 


, E (Lean, x 27; 


e domina estrictamente se 


Ef sn x; 


> Y; 


Se nemtx,,....x,) domina (y,, .... y,) nem (y,, ...Y,) 
domina (x,. .... x, ) então dizem-se dois elementos não 
comparáveis, 


A2.3.7 — Producto de Reticulados 


Sejam dados dois reticulados (O),,v.A) e 


(23, 4). Forma-se o producto cartesiano 2, X £)3. Con- 
sideremos um conjunto GC 2,X 23. Se quaisquer 2 
elementos de G forem comparáveis entre si, G é orde- 
nado pela relação de ordem dominância (estricta ou 
lata conforme o caso). 


A2.3.8 — Cadeias Maximais 


Empregando a relação de ordem w definido num 
reticulado é possível partindo de (0) ou (V) atingir 
respectivamente (V) ou (OQ). Escolhendo entre todas 
as cadelas assim construídas aquelas que não tenham 
clementos comuns (com excepção de O e da V) obtém- 
-se o conjunto de cadeias maximais do grafo. 


Esta definição autoriza que existam mais do que 
um conjunto distinto de cadeias maximais. 


A2.3.9 — Exemplo de um producto cartesiano de reti- 


culados 
G=V B Ev 
a, 
is Po: O 
Qo=0 


fla Ng 


Forme-se o producto cartesiano 2, X 25. Os 
pares que pertencem a esse producto são: 


(ago Bo) (o Bo) (as, Bo) (ara Bal 
lag B4) (a, B4) (as, B,) (args B,) 
(ao, Bo) (a, Bol (as, B.]) (as, Ba) 


Vamos usar o critério de dominância (=) —não 
estricta —, 


(ag, Ba) log By) > la, Bo) >, Bol>la, Bo) > lago Bo) 


O conjunto G € (2, X 27) formado por estes 
pares está ordenado pelo critério da dominância não 
estricta ( >). 


Note-se: 


i) Se 2, e 23 forem totalmente e estrictamente 
ordenadas, então obtém-se um reticulado vectorial 
(ver A2,4,3). 


ii) Se 2, e 23 forem distributivos também 


GC (2,xX“g) será distributivo 


A2Z2.4 — ESTUDO DE ALGUNS CONJUNTOS RETICULA- 
DOS MAIS USADOS 


24.1 — Conjunto de Zadeh ([0, 1], %,%]) 


Verifica-se facilmente que ([0,1],:,Y) tem as 
seguintes propriedades: 


TÉCNICA 438 


a) Fechado 

b) Associativo 

c) Comutativo 

d) Existência e unicidade de supremo e infimo 
e) Idempotente 

f) Abso:vente 

g)  Distributivo 

h)  Pseudo-complementado 


i) Aplicam-se os Teoremas de De Morgan 
« 5» significa «o mínimo de» 
cv» signífica o emáximo de» 


Neste conjunto [0,1], e para qualquer elemento 
a, & 2 =[0,1] é possível definir duas operações in- 
ternas fechadas: o producto e a soma disjuntiva. 


A2.4.1.1 — Producto (*) 
a, * q=aj=a . a 


onde (.) é o símbolo do producto vulgar. 
Esta operação interna (*) definida em [0,1] é: 


— fechada 

— associativa 

— comutativa 

— «O» é o elemento absorvente 


-— qi» é o elemento neutro 


A2.4.1.2 — Soma disjuntiva [(T) 


a;T a; = em e ta-—a nó 


onde 
«+» é a soma usual 
«a—» é a dife-ença usual 
e» é o producto usual 


Note-se que sendo a, e a; E [0, 11. 


; E E a &; 
donde o símbolo (--) pode aplicar-se e daí resulta 
dt E F0,1]. 


A soma disjuntiva é uma operação interna definida 
em [0,1]: 


(1) Ver exemplo em A2.3.9 
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— fechada 
— comutativa 
— associativa 


— O é o elemento neutro, 
AZ4.1.3 — As operações (*) e [T) são distributivas 
A2.4.1.4 — Tem interesse escrever a expressão de apli- 
cação sucessiva da soma disjuntiva: 


a TaT..Ta,= 


T u 
=(—-1) 3 a; +(— 1) pa la; e) 1. + 


+ (= qt la; “sas + E ) + 


au 
da 
a t=1 k 


l, a. 


+ 4-8] (e, .0,.... 04) 


A2.4.1.5 — Pode definir-se o operador (—), significando 


aj —a, à diferença vulgar, mas apenas se 


aj > q; 


A2.4.2 — Conjunto de Bool. 


Para este suporte 2 =(0,1) o reticulado ((0, 1), 
Rob) é: 


a — fechado 

b — associativo 

c— comutativo 

d — Existem supremo e Ínfimo, sendo únicos 
a — idempotente 

f — Absorventes 

g — Distributivo 


h — Complementado: 
à (0,1)=0 e v(O,1)=1 


i-—— Aplicam-se as leis de De Morgan 


AZ.4.3 — Reticulado Vectorial [!) 


Sejam dados n conjuntos A,, À.,.... À, qualquer 
deles completamente ordenado pela relação de ordem 
estrita (>) — a«>B significaa-f — e tendo um 
supremo e Ínfimo. 

Forme-sa o producto cartesiano 
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no 
XY <- X A, 
i=1 
Se usarmos a relação de dominação assim definida: 


(a, 08, ) >1b, o... Do) 
com lap... ap), Ibyosb,) E À 
se a, >b; Vie Et ssh 


então o conjunto 07% um reticulado distributivo mas 
não complementado. Com efeito: 


a) — Fechado 
ç El....a, PR PM Pa e Rj e Emote Pee) |] 


onde j 1€ =YV (a; ; D; |) e como A, é completa e 


estritamente ordenado (c, = a,) ou (c, =b;). 
Identicamente para (4). 
b — Associativo (*) 
c — Ceomutatividade (*) 
d— Todos os A, têm um supremo e um iÍnfimo 
e — Idempotentes (*) 
f — Absorventes (*) 


g — Distributividade 


g — Distributividade 


P[ (esse) A Uibysco), 


o 
ti 
c— 

[| 
| 


=A | v ti Apvady ucaby cva)lif 


TE nobre farei] [= 
= aTt..., Db, ia C; + ss] = (..., 9; utad 
em que 


g,=4lb,.c;), Vi 


v GA é PRP P sia 9; + «--1] —— RR hj...) 
hp =vla,.A lb; 1 
Por outro lado 


4 caspa ade fora Bio ed] = bica fito) 


=vla;.b ). Mi 


(*) E óbvio a partir da definição. 
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m,=v la; .c; ) Mi 
e 
io Bates (O Pe Mm; ER É — Damas hn, , csid 


n=4a(l.m )=Alvia,,b )vla,,e )] 


mas A, é completa e estritamente ordenado e portanto 
distributivo e dai 


h,=v [a,,4 (b,c,)] = 
=Alvla.b)vtap.ce,) =h; 
h — Não é complementado em geral a não ser que cada 
Ao seja (conjuntos de Bool). 
24.4 — Reticulado Lexicográfico 


Será um reticulado vectorial com os elementos 
totalmente ordenados. 


ANEXO 3 


A3.1 — Funções ordinais de um grafo vulgar sem cir- 
cuitos. 


Seja GC 46”, finito, anti-simétrico e sem circuitos. 


WY a família de todas as aplicações (funcio- 
nais ou não) de Y% em JL 


Y E Wum membro da família y. 


yl a aplicação inversa de y. 


Define-se, de um modo geral, N, |, um sub-con- 
junto vulgar, do seguinte modo: 


p—1 
| E *1 ç 
N; -| Uva (u;) E A Vc 
onde 


uewe O<k<p-1 Kkp En 


=| uva (u,)= O | 


e p&N, mais pequeno, tal que: 


Pa IN,)=0O 
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Os sub-conjuntos (vulgares) N, ,i & [0,...,r) for- 
mam uma partição de Ji e estão estrita e totalmente 


ordenados, isto é: 
(k<I) ES(Ng <N) 


Função ordinal 
definida como segue: 


Nivel de partição é dado pelo índice da partição. 
Assim, N, tem o nível ou ordem) k. 


Q (u,) =k se u, & No, 


A3.2, — Extensão da noção de função ordinal a grafos 
com circuitos 


Para o efeito formam-se as classes de equivalência 
do grafo G dado, onde a relação de equivalência é: 
«existe um caminho de u, para u e deu, pira u, ». 
Os conjuntos assim formados e máximos formam clas- 
ses de equivalência e constituem uma ordem total ou 
parcial. Se fo; total teremos a função ordinal desejada. 
Se for parcial construa-se o grafo vulgar sem circuitos 
a partir dessas classes e a função ordinal desse grafo 
sem circuitos será a função ordinal desejada. 

Na figura está representado uma parte da um 
grafo onde existe uma circulação entre os pontos u,, Us,, 
as. Us 


FIG. 1 


Esta circulação pode descrever-se pela equivalên- 
cia referida ao princípio: «existe um caminho de u, 
para u, e de u para u, ». 

Não sendo possível declarar quem precede quem, 
os quatro elementos u,, U,, U,, U, são tratados como 
se se «fundissem» num único elemento uZ e a figura 
acima pode ser substituída por esta outra 
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dum grafo vulgar sem circuitos é 


FIG. 2 


Desta forma suprimem-se os circuitos e então 
cai-se na situação geral de grafos sem circuitos. 
Por exemplo, na fig. 3 


1,2 8 3 é um circuito (1) 
9, 10 e 11 é um outro circuito (ll) 
4,5, 6 e 7 é outro circuito (Ill) 


7, 6 e B não é circuito 


FIG. 3 


Na fig. 4 está repreentado a forma simplicada e 
sem circuitos do grafo dado. 


FIG. 4 
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Education and research in Mechanics in the 


Federal Republic of Germany () 


RESUMO 


Considerações sobre os sistemas de ensino e in- 
vestigação da Engenharia e da Ciência na Alemanha 
Federal, a propósito do tema central da Mecânica, 
apresentadas na sessão sobre ensino e investigação do 
1.º Congresso Nacional de Mecênica Teórica e Aplicada 
(Dezembro de 1974). 


| am not quite sure of what was expected of me 
when | was Invited to talk at this Symposium. Therefore, 
vou just have to be satisfied with what | am going to say, 
even it |s not exactly what you would like to hear. 


| want to divide my talk into two sections: First 
| shall say a few words on education, and then | will 
present a very brief survey over research establishments 
that are active in the field of mechanics in West 
Ge:many. 


Let me start with a rather pessimíistic note. If you 
look back at the last 70 years of the history of me- 
chanics you will note that at the beginning of that 
period Germany was one of the undisputed leaders in 
education and research in mechanics. This was a time 
when mechanics was still deemed the indispensabla 
basis of engineering education, and it was a tima when 
many mathematicians still considered mechanics to be 
«the paradise of the mathematical sciences», as Leonardo 
da Vinci once said. 

Today, at the end of these 70 years, the picture 
has changed considerably: On the one side, mathema- 
ticians are ever more turning their backs to applications; 
they pursue fashionabla abstractions, unimpeded by 
tho messy realities of applied problems. On the other 
side, as a reaction to that development, there Is a 
growing tendency of our engineering departments to 
reduce the role of the theoretical sciences like mathe- 
matics and mechanics in the engineering curricula in 
favor of more «relevant» subjects, which are thought 
to be of more immediate applicability. Hence there Is 
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a growing gap between mathematics and engineering. 
This gap could, and should, be bridged by mechanics. 
However, tha desire to span the widening gap, puts 
more and more strain on mechanics. Some of my 
colleagues, who are still more pessimistic than | am, 
ara even afraid that mechanics might not survive, as 
an independent subject in the engineering curricula, a 
further increase of the gap. 

These remarks refer to my country, Yet, | think 
that they also describe the situation In some of our 
neighbor countries, To soma degree the tendency | have 
described seems to be universal. 


EDUCATION 


Let me now turn, after these melancholic remarks, 
to the education of peopla who are going to be active In 
mechanical science. By that | mean people working on 
research problems or research — related development 
problems in solid or fluid mechanics. Traditionally, a 
large part of these people has an engineering education 
as background; to a lesser, but still considerable degree, 
they were educated as mathematicians or physicists. 
Take as an example my professorial colleagues at the 
Institute of Mechanics of the Technical University of 
Darmstadt: We are eight professors; two of us were 
educated as mathematicians, two [including me) as 
physicists and four as engineers. The numbers at other 
technical universitias are similar. 

Because engineering is the main source for pro- 
viding research workers in mechanics, let me explain, 


(*) Paper presented at the ist Portuguese National Congress for Theoretical and Applied Mechanics. 
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at least cursorily, our system of educating engince:s. 
This might be of some interest to you, in particular 
since at the present time a change of that systems is 
under way. 


Traditionally, we had a two-lane road of enginee- 
ring education, The first lane leads to the academically 
educated «Diplom»-engineer, the second lane to the 
«Graduate» engineer whose education has a more pra- 
tical orientation. Young people who want to become 
«Diploms»-enginee:s have to change after four years of 
compulsory primary school to a secondary school 
(cObersculey or «Gymnasium»). After attending secon- 
dary school for nine years and after successfully passing 
the final examination of that school, the young people 
are allowed to study any subject at any German univer- 
sity. Particula:ly, they are entitled to study engineering 
at one of our technical universities or at any other 
university which has an engineering department. 


Here, | should explain that we used to have two 
types of universities in Germany: universities proper 
and technical universities. Both types of universitias 
have in common, among other things, that they are the 
only institutions of higher education which are entitled 
to award doctors degrees. The traditional universities 
teach law, economics, humanities, medicine, mathema- 
tics and natural sciences, but no engineering subjects. 
The traditional technical universities teach, at tha same 
academic level, natural sciences, mathematics and 
enginesring, but no humanities and no medicine. Inci- 
dentally, the position of mechanics between mathema- 
tics and engineering is reflected by the fact, that at least 
one traditional university, the University of Gottingen, 
has a chair of mechanics since the beginning of this 
century, though it has no other engineering chairs. 

Without going into details, | might mention that 
the difference in subjec matter covered by the two 
types of universities is gradually disappearing: A number 
of the traditional universities, and nearly all of the 
unive-sities founded after the war, are now taking up 
engineering. Likewise, the technical universities are 
increasingly turning to the humanities and to medicine. 

After 5-6 vyears at a technical university the 
student passes a final examination and takes the degrees 
of «Diplom»-engineer. If he wants so, he may stay for 
2-4 more years, doing graduate work with the aim of 
attaining a doctor's deg-ee in engineering. Most of the 
research workers in mechanics come out of this group 
of students. Tha education in mathematics and physics 
runs along quite similar lines and need not be explained 
here. 

Now, some words about the second lane of 
engineering education: Here the student usually first 
completes primary school at the age of 14-15. Then 
he learns, for 3-4 vears, a trade as an apprentice, at 
the same time attending a special trade school, After 
successful completion of trade school and apprenti- 
ceship, the young fellows go on to one of the many 
engineering schools, traditionally called «Technikum». 
At these schools they receive a sound training with 
much emphasis on the pratical aspects of enginsering 
and less stress on the theoretical side. Of course, a 
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deeper theoretical foundation of engineering can not 
be provided at these schools, because the students 
lack the training in mathematics that is (or was) supplied 
by our secondary schools, — After 3 years the student 
has to pass a final examination and then he leaves 
the engineering school as a «Graduate»-engineer. | 
should say that in the past a small number of them 
went on to a university, studying engineering there and 
finally receiving the degree of «Diplom»-engineer. 


Now, experience has shown, that Graduate-engi- 
neers are often doing the same type of work as Di- 
plom-enginesrs, for example in industry, but also in the 
civil service. However, they are oíten paid a smaller 
salary. And this is one of the more important reasons 
that trouble started in our system of engineering edu- 
cation. — | should say, that another reason was tha 
question of compatibility of engineering degrees in the 
Common Ma:;ket Countries; this is the same type of 
problem as already mentioned by Professor Zien- 
kiewicz. — Of course, it is only rational and just, that 
the same type of work should be paid with the same 
amount of money, Therefore, it would also seem ratio- 
nal and just to abandon a system of payment which is 
based on a formal qualification of tha employee, na- 
mely on his academic degree, and not on the type of 
wo:k he is doing. However, in my country, the gospel 
of tha equality of all men is inte:preted in a rather 
strange way, and we are told that the unsatisfactory 
situation can be changed only by giving everybody 
tha same type of education, regardless of his capabi- 
lities and, perhaps, even regardless of his desires. 


Às a consequence of this rather amazing theory, 
at the present time the official educational policy of 
the Federal Republic and its member states tries to 
merge the former engineering schools with ths univer- 
sities, thereby forming so-called «Comprehensive Schools 
of Higher Educations («Gesamthochschulenh). Six of 
these merged schools have already been established. 
As a matter of fact, this change goes parallel to a 
coriesponding change in the primary and secondary 
school system, but it would lead too far astray to 
explain this here. 


In these new comprehensive superuniversities all 
types of engineers will be educated, and all of them 
will receive the same titla of «Diplom»-engineer. Yet, 
funny as it may sound after what | have said about 
the reasons for that change, there will still be a diffe- 
rence between two types of «Diplom»-engineers: It is 
true, that during the first two years all engineering 
students will take tha ame courses. But then comes a 
separation: One part of the students will go on for 
only one further year, taking courses with a practical 
orientation. After a final examination they will become 
«Diplom-engineer of the first kind»; these correspond 
roughly to the former Graduate engineers. The other 
part of the students proceeds for two more years, 
taking courses with mora emphasis on tha scientific 
and theoretical aspects of engineering. They become 
«Diplom-engineers of the second kinds; these corres- 
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pond to the former Diplom-enginaeers. If they want so, 
they can complete their education with a doctor's degree 
after that. (*) 

Now, it is quite clear that this new system of 
engineering education has obvious and serious disadvan- 
tages. Here | want to mention only one drawback 
which has a special significance with respect to the 
education of people who are going to becoma activa 
in mechanics research: The coeducation of the two 
types of engineers mentioned during the first two years 
will necessarily be accompanied by a deterioration of 
the scientific and academic level of the courses in the 
basic sciences, mechanics included. This is so for two 
reasons: The first reason is the inhomogeneity of the 
student body. The student body will include a large 
number of students who have neither the ability nor 
the desire nor the need to study the more theoretical 
aspects of engineering. The contents of the courses, 
however, have to satisfy the needs of the majority of 
the students, Therefore, those students who want to 
proceed after the first two years with the more advanced 
courses leading to the «Diplom-engineer of the second 
kind» will do so with insufficient preparation. —- The 
second reason is the extreme inhomogeneity of the 
faculty. Besides the former research-motivated univer- 
sity professors the faculty will include the former en- 
gineering school teachers. These are usually engineers 
with some practical experience, but in most cases with 
neither the capability nor the readiness to do research 
work or any other academically advance work. It is 
quite clear that their standards will influence the level 
of the courses. It is interesting in this connection to 
note that the strongest opposition to the creation of 
comprehensive universities comes from the former en- 
gineering school teachers who are — justifiedly — 
afraid that they will be second class teachers only 
at the new universities. 

Without doubt, this development is a serious dan- 
ger for mechanics, because the continuation of mecha- 
nics as a science depends on a steady supply of 
academically educated engineers with a high quality 
grounding and interest in the theoretical aspects of 
engineering. 

They are many people today who think that 
«change» is synonymous with «progress». Yet, change 
does not necessarily mean a change to the better. | 
think that the creation of comprehensive universities 
is a change in the wrong direction. À system of edu- 
cation that neglects the truth that people differ in their 
abilities, desires and expectations cannot claim to be 
really progressive. 


RESEARCH 
The research establishments active in mechanics 


may be divided into three categories with different 
organizational structure: 


1. University institutes, most of them at technical 
universities. Support is provided mainly by the govern- 
ment of the state in which tha university is located, 
but a sizable part of the money for research is also 
provided by the German National Science Foundation 
(Deutsche Forschungsgemeinschaft, DFG) and a smaller 
part is earned through contract work, The German 
National Science Foundation is an independent funding 
agency for promotion and support of scientific research. 
lt derives the largest part of its income f:om the federal 
government and the state governments and a smaller 
part from industry. The Foundation distributes the 
money, after a careful review of research proposals by 
qualified outside scientists, with the purpose of suppor- 
ting specific research projects as well as supporting 
larger groups of scientists working jointly on interdisci- 
plinary projects, mainly at the universities, 


2. Industrial research institutes; many of them 
are part of big industrial companies, others are run 
jointly be several companies. Support comes nearly 
exclusively from industry. 


3. Independent research institutes. These are 
supported partly by the state governments and the 
federal government, partly by money from indus- 
try (through contract work) and partly by the National 
Science Foundation and other foundations. 


Besides these three types of institutes there are 
also governmental research establishments, However, as 
these play no prominent role in mechanics research, 
let me omit them. Also, because | am not really familar 
with industrial research organizations | shall disregard 
them. Let me only mention that traditionally the aircraft 
industry was at the forefront of research in mechanics, 
and this is still so. Of course, today, the term «aircraft 
industry» is meant to encompass the complete aero-space 
industrial complex. 

At the present time, there is a growing tendency 
to reduce the research activities at the universities and 
to favor am increased teaching load, Therefore, tha 
contribution of the technical universities to research in 
mechanics may be declining at the present time; but 
it still represents the largest and most important contri- 
bution. At every technical university, and also at those 
other universitios which have engineering departments, 
you will find a group of professors, together with their 
assistants and coworkers, who are responsible for tea- 
ching tha basic mechanics courses for all engineering 
students. These are also the people who do research 
in mechanies in the narrower sense. Usually from 5 - BO 
of these people are working together in special univer- 
sity institutes of mechanics. — In order to make clear 
my terminology, let me say that the universities are 
subdivided into departments, and that the departments 
are again subdivided into Institutes. This is a little sim- 
plification, but it suffices for the present purpose. 


(*) Note added during proofreading: Experience has shown that nearly all students of the comprehensive universitias try to get 
the higher qualification of Diplom-engineer of the second kind, because they think that this Increases their chances for higher pay 


after their graduation. 
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Traditionally the university institutes of mechanics 
ara heavily loaded on the side of solid mechanies: The 
research subjects range f;om rigid body dynamics and 
vibrations to elasticity, visco-and thermoelasticity, plas- 
ticity, fracture mechanics etc. Research might ba expea- 
rimental or theoretical or both. Only at a very small 
number of universities, including my own university, 
fluid mechanics is represented within these institutes. 
The institutes of mechanies may form an independent 
department of mechanics, as is the case at my univer- 
sity. At other universities different institutes of mecha- 
nics may be part of different engineering departments; 
this is perhaps the most frequent case. At still other 
universities the mechanics institutes may be part of the 
department of mathematical and physical sciences. This 
was the case at my university, before a reorganization 
split up the former big departments into a number of 
smaller ones. 

Due to the historical development, the teaching 
and research in fluid mechanics is usually done by sepa- 
rate institutes of fluid mechanics. Most of them are 
part of the mechanical engineering department. The 
dichotomy between solid and fluid dynamics is certainly 
not justified in view of the basic unity of mechanics 
and the didactic advantages to be derived from a unified 
approach to both solid and fluid mechanics. But despite 
all tha sweeping reforms going on at the present time, 
it is extremely difficult to change such anachronistic 
traditions! 

Of course, mechanics as a subject has no clearcut 
boundaries. So you will find in the engineering depart- 
ments many other institutes where research in mecha- 
nics in a wider senso is done. | mention the aerody- 
namics and aerospace institutes, the hydraulics institutes 
and the institutes of structural mechanics in the civil 
engineering departments. 

At the present time we have, in West-Germany, 
28 universities, 9 technical universitias, 6 comprehensive 
universities and 2 military academics with university 
status. Four universities and one comprehensive univer- 
sity are on the planning stage. At about 15 technical 
universities and other universities mechanics is repre- 
sented in some form. Not all the subdivisions of mecha- 
nics are represented in equal strengh at all the places 
indicated. As | am a fluid dynamicist myself, let ma 
point out that rather strong centers in fluid mechanics, 
including gas dynamics and rheology are to be found at 
the technical universities of Aachen, Berlin, Brauns- 
chweig, Darmstadt, Karlsruhe, Munchen, Stuttgart, and 
at the universities of Dortmund, Freiburg and Gôttingen. 

In the category of tha independent research Insti- 
tutes there are, apart from a few smalle; institutes, 3 
big research establishments. All of these pursue, among 
other topics, also research in mechanies: A. Tha Max- 
-Planck-Soclety, B. Tha German Aero-Space Ressarch 
Establishment (Deutsche Forschungs-und Versuchsans- 
talt fur Luft-und Raumfahrt, DFVLR), C. The Fraun- 
hofer-Society. 


A) The most famous of these organizations is 
tha Max-Planck-Society, which was founded in 1911 
under the original name «Kaiser-Wilhelm-Geselischaft». 
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At the present time the Max-Planck-Society comprises 
48 research institutes, nearly all of them in the field 
of natural sciences. These institutes are spread all over 
Ge:many. Of these 48 institutes only the Max-Planck- 
“Institute of Fluid Dynamics in Góttingen is relevant 
in the present context. It was established in 1924, and 
its first director from 1924 till 1946 was the famous 
Ludwig Prandtl. He was succeeded by Albert Betz, a 
wellknown asrodynamicist and longtime coworker and 
friend of Ludwig Prandtl. After Betz retired, Walter 
Tollmien became director; he is known first of all for 
his pioneering work on the stability of laminar boundary 
layers. After Tollmien's death, about 5 years ago, the 
instituto was reorganized. To the classical fluid dyna- 
mics and gas dynamics department a strong and active 
research group working in the field of molecular beams 
and low density flows was added. — It is certainly 
unnecessary to elaborate on the importance of the 
Max-Planck-Institute of Fluid Dynamics in the recent 
history of fluid mechanics! 

There is only one further Max-Planck-Instituta 
which | should mention here, because its subject area 
is at least adjacent to mechanics, and this is the Ins- 


titute of Plasma Physics in Munchen. 


B) The second big research organization is the 
German ÃAsero-Space Research Establishment (DFVLR). 
The DFVLR consists of 35 institutes plus 5 computing 
centers, spread over West-Germany. Contrary to the 
Max-Planck-Society the DFVLR includes many institutes 
which are activo in solid and fluid mechanics, Tha 
whole establishment is divided into & groups of ins- 
titutes: 


Group |: Fluid Mechanics, with the following six 
institutes: 1. Aerodynamics (Braunschweig); this ins- 
titute has a long tradition, going back to the twenties 
and being connected with names like Busemann and 
Schlichting. 2. Fluid Dynamics and 3. Dynamics of 
Rarfied Gases (Góttingen). These two institutes grew 
out of the former aerodynamic research institute in 
Góttingen. This aerodynamic instituts was originally 
founded by Prandtl who was also lts first director, later 
to be succeeded by Betz. The instituto was (and still 
is) located on the same site as tha Max-Planck-Institut 
of Fluid Dynamics. Both institutes worked always clo- 
sely together, some of the scientific staff, in particular 
Prandtl and Betz themselves, holding double appoin- 
tments in both institutes. 4, Applied Gasdynamics (near 
Kôln). 5. Applied Mathematics and Mechanics (Frei- 
burg): this instituta is headed by Henry Górtier, who 
is president of IUTAM at the present time. The subject 
area of the institute is mainly mathematical fluid dyna- 
mics. 6. Institute of Turbulence Research (Berlim). 


Group ll: Flight Mechanics, with the following thres 
institutes: 1. Flight Mechanics (Braunschweig). 2. Flight 
Control (Braunschweig). 3. Dynamics of Aero-Space- 
-Vehicles (Munchen), 


Group |ll: Properties and Strength of Materials and 
Structures, with the following Institutes: 1. Alrcraft 
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Construction (Braunschweig). 2. Aero-elasticity (Góttin- 
gen); this institute too is an outgrowth of Prandtl's 
original aerodynamic research institute. 3. Materials 
Research (near Kôln). 4. Methods of Construction 
(Stuttgart). 


Group IV: Propulsion and Energetics, contains seven 
institutes, for example: institutes of plasma-dynamics, 
reaction kinetics, rocket propulsion, lubrication. 


Group V: Electronics. This group contains three 
institutes, all in the Munchen area. They are working, 
among other things, cn control of and communication 
through satellites. These institutes cooperate with other 
European institutes and organizations responsibls for 
operating the Eu-opean and the German research satelli- 
tes, wich have been put into orbit with US rockets. 


Group Vl: Space Physics, Space Simulation and 
Space Medicine contains three institutes, 

The computing centers and a number of smaller 
institutes complete the DFVLR. The DFVLR has a total 
staff of about 4000 people. Among them is a lot of 
young and gifted and highly motivated scientists. Never- 
theless, at the present time he DFVLR has some trouble 
with regard to its future development. These troubles 
are more or less caused by the lack of a clear direction 
in German and European asro-space research. Fortuna- 
tely, this trouble does not affect too seriously those 
member institutes of the DFVLR which do fundamental 
research in mechanics, like the fluid mechanics iInstl- 
tute in Góttingen and Górtler's institute in Freiburg. 


C) The last of the big research establishments is 
Germany which | want to mention is the Fraunhofer- 
-Society. The Fraunhofer-Society is devoted to applied 
science, Todak, it consists of 17 institutes plus 10 
socalled «working parties», i.e. smaller institutes. | 
mention only the institutes in the broad field of me- 
chanics: 1. Solid mechanices (Freiburg); the institute 
works mainly in fracture mechanics. 2. Ernst-Mach- 


TÉCNICA 438 


“Institute for Shock-Wave Research (Freiburg). This 
institute works closely together with the Franco-German 
Research Institute at St. Louis (near tha meeting point 
of France, Germany and Switzerland). The Franco- 
-German Institute is a jointly government establishment 
with main research emphasis on high speed mechanics, 
ballistics, hypersonic flow, detonation and combustion. 
3. Strength of Materials (Darmstadt). There are wor- 
king parties for: Space research (Freiburg), applied 
materials research (Bremen), ballistics (near Basel), 
and others more. 

| am aware of the fact that | could present only a 
very incomplete and sketchy picture of the German 
mechanical research landscape. However, tha restriction 
of time forbids me to go deeper into tha subject. 

Perhaps | should concluda by saying a few words 
about organizations that support visits of foreign scien- 
tists and scholars to German universities or research 
institutes. Here | have to mention three organizations: 


1. The Alexander von Humboldt-Foundation, 2. 
Tha German National Science Foundation (DFG), 3. The 
German Academic Exchange Service (DAAD). 

The most important of these is tha Alexander von 
Humboldt-Foundation, with a tradition of more than 100 
years, At the present time the foundation supports visits 
of qualifiod young foreign scientists (between 25 and 
45), with complete university education, for periods 
of 1-2 years. The number of grants per year is about 
500. The grantees are paid a salary which is compa- 
rable with the average salary paid to German scientists 
of the same age group. 

The German National Science Foundation supports 
mainly guest professorships. The Academic Exchange 
Service supports qualifisd foreign students, below 32 
years of age, who want to study at a German university 
or research institute. Besides that, it supports visits of 
foreign scientists of less than three months duration. 
The embassies of my country provide detailed informa- 
tion on all the organizations mentioned and on tha 
proper procedure for application. 
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Continuous Hydrological Time-Series Discretization (*) 


ABSTRACT 


The available information about many continuous 
hydrological time series (XJ is often given by continuous 
time records with given finite duration and in many 
cases (i.e. rainfall and runoff records) the autocovariance 
function of X, can be approximated by a function or 
exponential type. Once the data have been deseasona 
lized, stationarity is usually assumed. A discretization 
technique is required to estimate the mean and the 
variance of X, which implies setting a value for the lag 
of time, L, between any two successive observations 
taken from the continuous record. The greater the time 
lag, the smaller is the number of used X, values, the 
weaker is the dependence between them and the shorter 
is the required computing time. The aim of this paper 
is to study the three effects to make possible the 
selection of an optimal value for L. Numerical results 
of general use are derived and presented. Obviousiy, 
these conclusions also can be applied to any nonhydro- 
logical iljme series providing the assumed statistical 
assumptions are acceptable. 


INTRODUCTION 


Considering the recent advances achieved in the 
ficld of recording systems for values of basic hydro- 
logic quantities (e. g., rainfall and runoff) it is expected 
that these records will ba a continuous series In the 
near future. These data are normally used to estimate 
tho time-series statistical parameters, among which 
mean and variance are the most important ones. Serial 
dependence is obviously a relevant feature in this 
estimation. 

The time dependence of hydrologic variables has 
been examined by many authors and the following 
general conclusions (3) can be drawn: (1) The auto- 
covariance function [R(/)] is usually significantly 
different from zero, and (2) once the data have been 


(*) Este trabalho fol 


Luis Valadares Tavares (*£) 


RESUMO 


A informação disponível sobre séries hidrológicas 
continuas, (X, |, resumem-se frequentemente a registos 
cronológicos com duração finita. Em muitos casos (por 
exemplo, para precipitações e escoamentos) a autoco- 
variância de (X, ) tem um andamento exponencial e 
a hipótese de estacionariedade é aceitável para a série 
dada desde que se tenham filtrado as suas compo- 
nentes sazonais, 


Para estimar a média e a variância de (X, ) é 
necessário discretizar os dados e, portanto, adoptar 
um valor para o intervalo de tempo, L, que medeia 
entre observações consecutivas lidas no registo conhe- 
cido. Quanto maior L, tanto menor é o número de 
elementos da emostra assim obtida mas tanto menor 
é a sua interdependência bem como o tempo e o tra- 
balho necessários ao seu cálculo. O objectivo deste 
artigo consiste em estudar estes efeitos de modo a 
poder optimizar L apresentando-se elguns resultados de 
uso geral. 

Como é óbvio, as conclusões obtidas são também 
aplicáveis a séries não hidrológicas desde que satisfa- 
cam es hipóteses admitidas. 


deseasonalized, stationarity sometimes can be assumed. 
Often, the type of the function, Ali), is one of those 
represented in Fig. 1. 


In Ref. 1, several processes with R(i) of types a 
and b are described. Some common examples are men- 
tioned In Table 1. 


Thus, an important problem in hydrology is to 
optimize the method of estimating the mean and va- 
riance of continuous time series with the autocovariance 
function of type a or b using continuous samples of 
given finite duration. To estimate these parameters, data 
are discretized after setting a value for he time lag 
(L) between any two successiva observations (X,, X, | 
taken from the continuous record. Thus, the selection 
of L must be studied. 


ublicado no JOURNAL OF THE HYDRAULICS DIVISION da AMERICAN SOCIETY OF CIVIL ENGINEERS, 


em Janeiro da 1975 se 49.63) e baseia-se no capítulo Il (Princípios e potencialidades da inferência Estatística) da disser- 
m 


tação «Métodos de 


ulação Hidrológica» que o autor apresentou para obtenção do grau de doutor em Ciências de Enge- 
nharia (Hidráulica) pela Universidade Técnica de Lisboa. 


(*+*) Centro de Sistemas Urbanos e Regionais das Universidades de Lisboa (INIC), IST, Av. Rovisco Pals, Lisboa. 
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TABLE 1 ——- Models for Types a and b 


MODEL TYPE OF R (1) 

(1) (2) 
Xe —un=S(X ,—ul+te; with & >0 a 
(X6— u) = 6, (X, 4 =) +&(X,o-u)+e: withb, >0 and &,>0 a 
4 — )—PMX y—u)= sm 0,8: with ?, >0,69,>0,and 6,/b, <1 a 
X-un=% (Xey-ul+te; with&, <O0 | b 
M-0)=8 (My +e XMey—u te; withg, <O0 ande, >0 b 
(Xe— a) — &; (X-g-un)=e—0,e, ; with Pb, <O and 6,>0 b 


As L is increased, the serial dependence between 
the discrete values is decreased so that each used 
values (X, ) contains a higher volume of information 
that is not given by other used values. Increasing L 
implies getting a smaller number of discrete values 
(N) if a finite duration continuous historical record 
is available, Obviously, this method requires a smaller 
amount of computational work. 

Therefore, the optimal value of L is the result of a 
compromise between the volume of information (1) 
extracted from the record and the computing effort 
spent to process that information (E) .Thus, | as function 
of L must be studied and is examined subsequently. 

The valua of E is afunction of N and its expression 
depends greatly on the working facilities and objectives 
of each particular case, Therefore, postulating a general 


FIG. 1 — Types a and b of Correlation Function [ P (i) |] 


expression for the function, E, is not very realístic. It is 
preferabla to analyze the sensitivity of / to variations in 
L. Then, for any particular case, the optimal range of 
L values can be determined easily. 

Many ways of measuring | have been proposed (2). 
Again, the most appropriated methods depend on the 
objectives of the study. However, all are functions of 
the variance (or the standard error) of the estimated 
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parameters, as it is easy to calculate their values if the 
variances are know. Thus, the variance is basic to the 
definition of ! and herein it will be studied as a measure 
of |. 


UNBIASED ESTIMATOR OF MEAN (u) AND ITS VA- 
RIANCE (V,) AS FUNCTION OF L 


lt is well known that 


N (1) 

E No 

such that E = and N= -— 
lu) =p L (2) 


in which N, = an arbitrary integer constant. 
To calculate tha VAR (u): 


V, =VAR (4) (3) 


| 1 N 2 
n=E|43 lts 8 | (4) 


Then 
1 N N 
he É B M-u+t H X—4) oq | 
N “51 LE) 
(5) 
1 N 
n= | Mete O mn 040 | (6) 
Let S, be defined as 
N 
= Y nN-j ] 
TYPE a—lf R (i) is of type a: 
p (i) = RI] (8) 
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with A = e 2”, By definition, L is equal to! if N=N, 
and a« = an appropriate constant, Thus 


” ade NR (1-RYIR 
SN Nº |141-R (1-R)? (10) 
or V, = Bo? (11) 
with A= a o BR Re 8 =RO)R 
“ NIS-R (1-—R)? (12) 
1+A 
= ———— 
N (13) 


lt is important to calculate the limit of V, as L 
tends to zero (LV,). Then 


LV, = lim V, (14) 
EQ 
1+A 
LV, = e? lim 
R Loo N (15) 
Because R = Tl = 1-al (16) 
1.0 
B 
E 
0.9 
0.1 0.8 
— (1.7 
0.6 
0.5 
0.4 
0.3 
0.2 
di 0.1 
type a 
0,001 ' 
| 2 20 200 N 
ij ||| — 
1 102 10) 19º 10º 
so 5 0.5 L 
10 1 0.1 0.01 0.001 


FIG. 3—B as Function of N (or L) and of a“? (type a) 
Assuming N, = 100 


one obtains 


2 L(1 —- al) 
LV, -| N, E 
A | dd o | e 
VS Eoo (1-—R)? (17) 
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ZaN,-—2 
LV, =——— q* 
' a? Nú (18) 
2 
or LV, = 
E (19) 


Tha values of B are represented in Fig. 2, assu- 
ming that 
N, = 100 (20 a) 


e7= 0.1; 0.2; 0,3; 0.4; 0.5; ...; 0.9 
(20 b) 


L = 10; 107; 101; 0.2; 0,3; ...; 0.9; 


1.0; 2.0; 5.0; 10.0: 20.0; 50.0 (20 Cc) 
TYPE b-—lf R (/) is of type b: 
1=0?B (21) 
with 
1+A 
RR 
N (22) 
and 
o 
Ã=—— 5 
No! (23) 
such that 
N—1 
S, = (9) AN — 1) | 
a" (24) 


Using L<.1 has no meaning any more because the 
correlation function of type b is not defined for L < 1. 
Thus, only L values greater than or equal to one will 
ba considered, 

An analytical expression for S, and S, also can 
be deduced for this type of p (/). By definition: 


s= 3 N—n(—- ey Sê 
N—1 sp 
S, = p3, IN — 1) 67% (26) 
Therefore 
NR? R? 
Sd Jp (-AR (27) 
and if L is even: 
NR R 
“TOA MAR (28) 
If L is odd: 
S=(N—-WI-R)+HIN—2)R?+ 
+IN—3)(— R3) +. + (— R)SSI 
(29) 
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Therefora 


N(=-R+HREERIE-R)+RAR— 2R? 


S, = 

; = (30) 
because R' =0 (31) 
assuming N,= 100; e *=0.1; 0.2; ...; 0.9; and 


L=1,2,3,..., 9, 10, 20. The values of B are p:e- 
conted in Fig. 3. 


UNBIASED ESTIMATOR FOR VARIANCE (9?) AND ITS 
VARIANCE (V,) AS FUNCTION OF L 


The value, o?, Is defined considering: 


TAM = N=-1 
A MN NE 


2 q? e | 
cs E Us 
a > NI) R | (32) 


Then 


N 
1Y a Did — 
| 3 Xj —a) | W>T-AdA am 


1.0 
05 
B 
0.1 
0.05 
0.91 
0.008 
type b 
2,061 Ep 
0.0005 = 
N 
a OR + 
5 19 20 40 so dO tm 
0. A 4 
ER Tá de IR A AR SR DR 
2 í 16 5 1 


FIG 3--B es Function of N (or L) and of e” (types b) Assu- 
ming N, = 100 
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Ds 
a - (34) 
º N-1-A 


in which A = the loss of degrees-of-freedom due to the 
serial correlation if tha variance is estimatad. 


To calculate V,, q = u is assumed. Thus 


V=———— vAR| W (X.- já 
PANA BU e ni 
== E VARIK,— qt + 
à = iÃ didi 
N 
+ >) COVAR [(X, — u)?, 0600) 


Therefore, a relationship between COVAR [X, — u?), 
IX; -— uj2] and COVAR [X, — 1), (X; — u)] must be 
derived. If X and Y are two reduced normal va- 
riables: 


03 
B 
81 
005 
8.01 
o.00s 
É 
0.001 Nhogo 
o 
ER! 
0.0005 1 0.90 
L 
kgs 
0.95 
No. 
+— | | — = 
5 18 Jó e] Eb BO 100 


FIG. 3 — Continued 


COVAR (X2, Y2) 


X2, Y2) = 
pi Adi Tea 2 (37) 
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- 3 = + Ee el y2m= + [a 
COVAR (X2, Y2) -| , f : (x — FS — 1) PF UXS, » Va) d (Xº) d (73) (38) 
+ — q 


COVAR (X2, Y2) Ta É sea e 
X2, Y2) = Xº — 4 FP = 1) Fo dt(y: FX) d (Xº 
f 3=0 [f 92=0 | Es | PRA (39) 


such that f (X2, Y2), f(Y2/X2), and f(X2) = the probability density function (pdf) of (X2 N Y2), Y? given X2, 
and Xº?, respectively, Then 


j x? = 
COVAR (X2, Fº) = f (X2 — 1) Lo(lX, Y)]2 (X2— 1) f (X2) d (X2) (40) 
7 3 
or 1 
| | O = +— mms 
COVAR (X2, Y2) = 2 [o(X, W)]? (41) N-A (46) 
Thus 
+ 28, 
p (X2, V2) = p2 4X, Y) (42) O ld AM (47) 
1.0 TYPE a—lf R (li) is of the type a: 
C 
0 type à pli)=R' (48) 
NR? (1—-R2N)Rº? 
G=[—— = 
0.01 2 1-R2 (1— Rº): (49) 
0.001 
-u 0.12 
0.00 B = 
0001 D 
0.10 
0.00001 0.1 
0.9 0.00 
0.000001 
2 20 200 N 
1 (A+ ———+————+— 
10 192 103 10* 10º 0.06 
50 5 0.5 E 
10 : 0.1 001 0.001 0.046 


FIG. 4— C as Function of N (or L) and of 0”? (type a) Assu- 


ming N, = 100 ud 
If X and Y are not normally distributed, but f(X, W) is 0.00 
approximately normal Eqg. 42 is still acceptable. Other- isa 
wise, Monte-Carlo simulation must be used to estimate D. 
the relation between p (X2, Y2) and p (X, Y). Assuming aos 
Eq. 42: 0.006 
2 os 2 0.004 
Fo —e —— ) q4 S 0.002 
NA Cas CO quai ps 
0.000 
with 
N—1 
S= DeX. Xp (N-i) | | Fra 
[1 (44) 10 01º 001 0.001 
or FIG. 5— C as Function of N (or L) and of e'Z (type b) Assu- 
1 ming N, = 100 
4 ese sda 45) 
In Eg. 49, Ram is negligible compared with tha other 
with terms. Thus 
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Rº 
Nat = ————— 
1—R? 
Penão 4 (50) 
And if 
L<<Nç(1—R?) (51) 
then 


2 Rº 
- 14204 —— 
TES (o e Er A (52) 


such that 4 is given by Eq. 12. 
Using N,= 100, Eq. 51 is acceptable for 
O14<R<OSIfL<<19(orND>D>õêõB).In this case 


N—A (53) 


FIG. 6-= C as Function of N (or L) and of e-z (type b) Assu- 
ming N, = 100. 
2 Rº 
De —— ea 
N— AJ — Rº) (54) 


The limite of V, if L tends to zero is also an 
important result. 


1 1 
LV? = lim Vo = 44 ( Rs, 


L+0 cima DR Ti dedo 
because 
Â 2 
Loo Nº an, (56) 
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end 


Ban: = = 0 (57) 


Therefore 
ZaN,—A 201 


=)93d-——————— e ——— 
Mi So E ai, (58) 


In Figs. 4 and 5 the values of C and D are presented 


fo: this type. The values of e “ and L are those assumed 
in Fig. 2. The value of D was calculated using Eq. 50. 


TYPE b— If R (i) is of type b: 


A] 


0.81 


E,001 


6.090) 


FIG. 7-—D as Function of N (or L) and of ez (type b) Assu- 
ming N, = 100. 


1 y 
E Vo = Cot+Da, (59) 
with 
1 
C=— 
N-A (60) 
28, 
“—qN-AR (61) 
2 
ia (62) 
N=—1 
“= ba, pt AN — i) 65) 
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N 
= ae — | 
S, pa) 2X, Xi) UN — dd) (eai 


in wich S, and S, were calculated using Eq. 28 or Egs. 
27 and 29, The C and D values are presented in Figs. 
6 and 7 using again NV, = 100; ez = 0.1; 0.2; ...; 0.9; 
and L=1,2,3,... 9, 10, 20. 


CONCLUSIONS 
For continuous stationary hydrologic tima serles 
with type a or b correlation functions the following 


conclusions can be derived; 


1. The unbiased estimator of qu is in any case the 
classical function 


N (65) 


2. To estimate the variance (o*), the serial de- 
pendence is responsible for an extra loss of degrees-of- 
-freodom (A): 


N ma 
2 (X, Es a)? 
um j=1 
E e aie—esmiãa 
N = — À (66) 


A can be calculated easily using Eqs. 12 and 23. 


3. For any particular time series the variances 
of u and o? are functions of the discretization process 
definsd by the L value: 


V, = VAR (u) = 02 B(L) (67) 
Vo = VAR (0º)=2 0! [CIL) + OD (L)] (68) 


This last formula implies X distribution being normal 
or quasinormal. 

The functions, B (L), CL), are presented in the 
Figs. 2-7. These functions are far from linear with L. 
For tha type a correlation function and introducing small 
errors, they can be approximated by 


100 
B(L) = -— In — * "HH L t, 
(L) = Pp — pn — nSbdt 
B IL) =p f L as Lp (69 b) 
L 
C(L)= ums A H sr Ê t 
(L)=pc-— q 100 c<L<te (70 a) 
100 
D(L) = Qntn— :s Hs L<t 
(L)=pp +7p L p“Sl<tp (71 8) 
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pb=rs;dté > Lá (71 b) 


in which Pp Jp Fo Sp tp Lp. Pre Ace fesSc te 
[] Le. Pp Tp Fon» Sp» tn» and Ln = positive para- 
meters the values of which depend on the value of 
e” *, Therefore, for any particular record: 


dV, o? 
a = O L t 
o 4,7 see Sê oa vi 
dv, 
—-— = O Ff LÊ Ê 
dL > Lp (72 b) 
dV, 20t 201 
dE SE CID 
if max (sc.sp) <L<min testo) (73 a) 
dv, 
—s = (je if É max (L-,lL 
dL > Leto] (73 b) 


Thus, the reduction rate of V, by decreasing L 
decreases with the inverse of L, and for L higher than 
L p + increasing N is quite irrelevant to decrease V,. The 
variation of V, with L depends on the values of « and of 
ot. In any case, for L > max (Lc. Lp), increasing N Is 
also irrelevant to decrease V.,. 

If the cost processing the data (E) is proportio- 
nal to N: 


dV, o? , 

“ge Proportional to — q» N | H sp<L<ta, 
(74 a) 

dv, 

-—— proportional to O; MH LS5L 

dar” 2 Lp (74 b) 

dV, ii 204 20+ 

es rt to — seg asian o 

a; Proporional to = dg N + qp N if 

max (sc.Sp)<L and L>minttoto) (75 a) 

“JE proportional to O; if L>max(Lç Ls) (75 b) 


Thus the effectiveness of increasing N to reduce V, 
or V, is either inversely proportional to N (Egs. 74 a 
and 75a) or equal to zero (Eqgs. 74b and 75 b). 


4. The limits of V, and V, as L tends to zero can 
be deduced analytically if 


e Ui) = RI (76) 
obtaining: 
' 20? 
LV, =— 
a (77) 
20t 
LV, = 
aN, (78) 
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in which M, = an arbitrary value of N for which L = 1 
by definition and « is defined as 


e=-InaR, (79) 


such that R, = the correlation coefficient between su- 


ccessivo X values, using N=N,. H Nº =KN, and 
K = any real constant, the definition of X implies: 


Ra= a (80) 
If tha correlation function is Eq. 76, then 
". pllk | 
Ro = Ru (81) 
and 
a No = q MN, (82) 


Therefore, (a N,) does not depend on the N, value. 


5. To use the presentei graphs a rough estimate 
of R (i) is sufficient, In hydrology, this estimation can 
be obtained fÍrom previous hydrologic studies, from 
general geophysical considerations, or from a pilot ana- 
lysis of recorded subsamples. 


6. Nowadays, important hydrologic magnitudes 
(e. 9. rainfail depth and runoff) are often recorded 
using continuous measuring systems. For each basin, 
these records ore known for a tima period with a given 
finite duration. Therefore, in hydrology, to estimate the 
mean and variance of a stationary continuous sampla 
with a given finite duration is a common and fundamen- 
tal problem. Often, the fractions of project work and cost 
that are required to carry out this estimation process 
are high because it is used repeatedly. Thus, impro- 
ving this process has important scientific and economic 
effects, To solve this estimation problem the continuous 
sampls must be discretized. Therefore, the decision 
variable is L, defined as the adopted discretization time 
interval, Then, improving this estimation process is 
achieved by studying the method of selecting L. Se- 
lecting L is making a compromisse between the esti- 
mation precision level and the required amount of work. 
Optimizing this decision implies obtaining the trade-off 
between measures of these two quantities. This rela- 
tionship was derived herein and is presented in Figs, 
1-7. These results can be applied easily to pratical 
hydrologic problems to avoid undesired precision levels 
as well as irrelevant increments of the estimation work, 


ACKNOWLEDGMENT 

This study was sponsored by the Calouste Gul- 
benkian Foundation through a research project directed 
by E. R. Arantes e Oliveira. 
APPENDIX | — REFERENCES 
[1] Box and Jenkins, Time-Series Analysis: Forecasting 


and Control, Holden-Day, San Francisco, Calif., 
1971. 


324 


[2] Fisher, R. A. The Design of Experiments, Oliver 


and Boyd, 1966. 


[3] Yevjevich, V., Stochastic Processes in Hydrology, 
Water Resources Publications, Fort Collins, Colo., 


1972. 


APPENDIX || — NOTATION 


The following symbols are used in this paper: 


1 


A,B,C,D,E 
COVAR (a, b) = 
E 1) 
to) 


j= 
ft = 

KlLLoilc Lo=- 

LV LV, 

au = 

No = 
PoPç Pp 76 7C7p- 
R,R* = 

Ro = 

Rs = 


Ri) = 


Fporerp= 
S4 92 5h: Sc: S p= 
q PAD PD 
Voo = 

VAR () = 

de= 


pla,b) = 
p Ui) = 


gé = 


= probality density 


real variables; 
covariance of a and b; 
expected value; 

fun- 
ction; 

volums of information; 
positive integer variable; 
positive real variables; 


= real variables; 


positive integer variable; 
particular value of |V; 
positive real parameters; 
positive real variables; 
particular value of R; 
particular value of R*; 


valus of autocovariance 
function of X for lag = i; 


positivo real parameters; 
positive real parameters; 
instants of time; 

real variables; 

variance; 


real variable; 


= value of X at instant t; 
= positive real variables; 


= random variable for ins- 


tant ft; 


moving average real pa- 
rameter; 


mean of X; 


= estimator of q; 


correlation coefficient be- 
tween a and b; 


value of correlation fun- 
ction of X for lag = i; 


variance of X; 


estimator of c?; and 


regressivo real 
ters. 


parame- 
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